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CHAPTER 1 


INTRODUCTION 

There has been a continuing interest in the use of silicon-on- 
sapphire, SOS, as a material for fabrication of high frequency electronic 
devices. The SOS technology allows dielectric isolation of devices by 
etching out isolated islands of silicon film on the high resistivity 
sapphire substrate In which individual devices are constructed by dif- 
fusion, oxidation, metallization, etc., implemented using photolitho- 
graphic definition techniques. Because of the low lifetime in thin 
silicon films on sapphire, practical bipolar transistors have not been 
realised using this technology; however, MOS circuitry Is now manu- 
factured in SOS (Inselek, RCA, and perhaps others.) It seems likely that 
CCD's (Charge Coupled Devices) built in SOS may be attractive for too 
reasons. ■ 

First, it may at some time be desirable to Integrate CCDs with other 
high speed SOS logic circuitry. As we shall see, CCD's will probably 
require thicker films, as will probably be the case for bipolar tran- 
sistors, and this poses technological problems for etching. The thin film 
is more desirable since the etch time is reduced and the photoresist 
protection problem is less severe. 

Perhaps a better case can be made for the desirability of fabri- 
cating CCD imaging devices in SOS. There are problems in accomplishing 
this, but the problems of building other CCD imagers must also be con- 
sidered In comparison* The central problem in fabricating a CCD is to pre- 
vent potential barriers (or wells) between the transfer electrodes. There are 
three techniques which have been successful in overcoming the problem of 


potential barriers in the gap which result in the case of a simple 
structure with gaps wider than 1-2 u m. (1) The technique used by Texas 
Instruments 1 makes the gap very small by depositing half of the aluminum 
electrodes and then anodizing them, producing an A1,,0 3 insulator, roughly 
1000 A thick, to insulate these electrodes from a subsequently deposited 
aluminum layer. (2) A technique proposed at Bell Labs 2 utilizes intergap 
doping to establish a potential within the gap to a value between those 
under the transmitting and receiving transfer gates. (3) Researchers at 
RCA 3 have published a number of paners giving results obtained with tXv'O 
levels of metalization, where one conductor level is doped polysilicon. 
Another technique which has been mentioned is the resistive sea of poly- 
silicon in which the transfer electrode pattern is defined by doping. A 
short discussion of these techniques will be beneficial. 

The first technique appears to be the best for solving the problem. 

It is simple in concept and implementation and produces, insofar as 
reported, no undesirable side effects such as oxide contamination for 
example. A high conductivity transfer electrode structure is defined in 
this way and the gap potential is fixed by the surface potentials under 
adjacent electrodes. Prom the standpoint of imaging, the channel is closed 
to incident radiation so that backside illumination is required. Therefore, 
the wafer must be thinned in the channel region to a thickness of roughly 
30 u n, i.e. , the spacing between adjacent cells. The second technique 
essentially results in a self-aligned bucket-brigade structure which hah a 
higher transfer inefficiency and consequently a lower upper cutoff fre- 
quency than a true charge coupled device. Otherwise, the structure allows 
wide gaps for enhancing the photon receiving area of the channel. 
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The third technique utilises doped polysllicon electrodes which at 
best have a resistivity over ten times that of aluminum. Consequently, in 
long structures the resultant RC transmission structure for the polysilicon 
transfer electrode structure produces phase delay and distortion which 
limits the upper frequency of operation. The polysllicon-oxide-silicon 
thin film structure will result in spectral selectivity of the incident 
radiation, and there will of course be some absorption of radiation in the 
polysilicon. This latter feature of absorption will be inherent in any 
scheme except number two above. 

Although all of these schemes may prove satisfactory in some appli- 
cations, they obviously have problems associated with them. In a preceding 
report, 4 it was suggested that a thin film silicon-on-sapphire structure 
might be feasible. From an optical standpoint such a structure would be 
superior because the sapphire is transparent and most of the photons would 
be absorbed in the depletion region of the CCD channel, Kranzer 5 has 
reported lifetime profiles in SOS films which show that the thermal 
(leakage) current from depleted regions near the silicon-sapphire interface 
would be excessive for sensitive SOS imaging devices. His results suggest 
that thicker silicon films will be required for successful applicat* on of 
SOS to CCD imaging devices. Either of the techniques (1) or (3) above 
should be suitable for transfer electrode structure, 

Kranzer 6 has done extensive experimentation with gate controlled 
diodes in which the capacitance and leakage current were measured for various 
conditions. Meyer 7 has also reported results from experiments with this 
type of structure. These results should be valuable for the design of CCD 
imagers and perhaps for bipolar transistors fabricated in SOS. 
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Figures 1.1 and 1.2 show typical results obtained by Kranzer. 6 His 
published data 5 was obtained using the one-dimensional Poisson's analysis 
to reduce both capacitance and diode reverse current vs. gate bias data to 
obtain lifetime and doping profiles for depletion depths short of the 
silicon-sapphire interface. The results of a one-dimensional analysis can 
be expressed as : 


X d " f si /C si 

1.1 

1 = 1/2 qHjA/OI r /aX d ) 

1.2 

W - - 3 V ax d 

1.3 

C sl - C T C ox /(C o* - C T> 

1.4 


Where C^, and C qx are respectively, the measured gate less 

overlap capacitance, the series depletion capacitance, and the gate to 

depletion layer capacitance through the oxide. Other variables are t , the 

carrier lifetime, n, the intrinsic carrier concentration, e . the silicon 

i si 

permittivity * V^, the gate voltage, I r , the diode reverse current, q, the 
electronic charge, A, the area of capacitance C and X^, the depletion 
layer depth. 

There are certain features of Kranzer's data 6 which are not explained 
by a one-^dimensional analysis. Furthermore, it would be of interest to 
have the results of a more complete analysis to study the data in the region 
in which the one-dimensional interpretation has reasonable validity. This 
study was conducted to obtain an. analysis based upon an essentially two- 
dimensional model. The structure studied has cylindrical symmetry so that 
only the Z-axis and radial coordinates are required. 
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Chapter 2 presents the one-dimensional analysis of the MOS structure 
as it is applied in the succeeding work. Chapter 3 presents a somewhat 
simplified two-dimensional analysis which is useful for interpreting the 
experimental results. Chapter 4 presents a more complete two-dimensional 
analysis. Chapter 5 gives the conclusions of this study. Appendices A, 

B, and C give the computer programs developed to implement the analysis. 


CHAPTER 2 


ONE-DIMENSIONAL DEPLETION 
LAYER ANALYSIS OF MIS STRUCTURE 

A diagram of the structure to be modeled is shown in Figure 2.1. 

The capacitance structure is annular with typical dimensions given in the 
figure where the aspect ratio of the differential radius, , to the 
thickness is greater than 100. If the silicon-sapphire interface is 
accumulated with mobile carriers, thus providing a conductive plane, then 
clearly the silicon depletion field variation across the film thickness 
can be closely approximated by a one-dimensional Poisson equation solution. 
The presence of the Pi inner annulus and the N+ outer annulus will give 
rise to a small radial field component. The fact that a rectifying barrier 
occurs at the junction of the P+ and N annulae will make the radial field 
component larger in the junction region and smaller elsewhere. Therefore, 
until the depletion extends to the silicon-sapphire interface, the one- 
dimensional approximation should be valid. 

Figure 2.2 shows sketches of the assumed charge, electric displace- 
ment, and potential distributions. Surface charges at the Si-Si0 2 and 
Si-Sapphire interfaces are assumed. The existence of the former is well 
established, and the latter is required to prevent depletion of the inter- 
face when a two-dimensional model is used. The two-dimensional calculation 
discussed in Chapter 4 showed that the field due to the junction would 
penetrate the sapphire and deplete the more lightly doped interface. If 
one assumes the thermally generated carriers in this region are collected, 
this gives rise to a leakage current, due to generation in this low life- 
time region, which is much in excess of that observed with small gate 
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voltages. A fixed positively charged layer will keep the interface in 
accumulation until the field component due to the gate electrode depletes 
this region. 

Application of the Poisson analysis to this problem is relatively 
straightforward. It is advantageous, however, to consider two cases, 
i.e., Case I, depletion is short of the Si-Saphire interface, and, Case II 
depletion extends through the silicon film. 

Case I. Depletion thickness less than film thickness. 
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Then the potential Is: 
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+ V* - W /E S i - *** 


For any depletion depth = x - t Qx , 


V G = " CQ B + Q ss2 )/C ox " V C si + 
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Equation (11) may be used to calculate the gate voltage required 
to deplete any arbitrary doping profile to the depth x^, Note also that 
from (11) : 


= - (1/C_ + l/C si )dO B + d(A^) 


= - dQ B /C T + d(Ai^) 


where is the measured gate capacitance. Carrying out (carefully) 

the indicated differentiations with respect to x^, one obtains: 

C T dV G 

N B <X d J " ” q d57 2 * ; 


Further, from (6): 


£ Bi U - V 


e. (1 - C„) 
x. = t 

d C ox C N “ 


where C „ - C T /C ox is the normalized capacitance which could be obtained 
from measurements. Note that equations (13) and (14) provide a parametric 
pair which may be used to estimate the doping prof •» la N g (x) for the 
depletion region giving rise to the capacitance C si which is included 
in the measurement. 

The leakage current due to generation of carriers in the depletion 
region of the field induced junction is simply given by: 


61 = qA 

gen ’ 


t +x , 
ox d 


2t (x) 


* less overlap capacitance. 
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2.16 


T<x d ) = |<m 1 /C3I gen /3 Xd ) 

Again (14) and (16) form a parametric pair for establishing the lifetime 
profile from junction current vs. gate voltage measurements if one assumes 
that the measured junction current, 1^, is equal to the generated current, 
I *^6 ure shows typical doping and lifetime profiles obtained by 

applying these parametric equations to TCranzer's data. 

Case II. Full depletion of the silicon film. 

In this case, only the regions for solution of (1-a) and (2-a) are 
extended. On the sapphire side of the si-saph. interface we have: 


D 

sa 

= D 

ox 

+ %sl + ’b + ‘W 

1-C 

sa 

= ifi 

xo 

- (D ox + <W /C S i - A * 

2-c 


When this full depletion occurs, the field generally will be determined 
by boundary conditions and geometry so that iJj x^ill no longer be a 
constant value across the si-saph. interface. However, (1-c) and (2-c) 
may be used to establish the interface boundary conditions for a field 
solution in the sapphire region. Combining (1-c) and (2-c) , we obtain: 

*b. ■- % + ^Bl - D sa> /C T + W C ox + V G - 2 - 17 

In the following, (2-17) is used in obtaining a numerical solution for the 
field in the sapphire when the through depletion occurs. 

The equations obtained here can be used two ways. As pointed out, 
the pairs (13) - (14) and (14) - (16) can be used for analysis of experi- 
mental data as was done by Kranzer. If the profiles t (x) and Ng(x) are 
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available, then (6), (7), (11) and (13) can be integrated (we use 

numerical methods) to obtain the capacitance voltage C T vs. (or 

C M vs, V„) and 61 vs. V n curves. 

N G gen G 


CHAPTER 3 


SIMPLIFIED TWO-DIMENSIONAL ANALYSIS 

This work was done after that repotted in Chapter 4; consequently, 
some insight into the electrical behavior of the QCD on SOS had been 
obtained. It was known that the 1-D depletion model gave a fairly 
good approximation before through-depletion if Q ss: j (the fixed 
charge at the Silicon-Sapphire interface, SSI) was large enough to 
prevent depletion from the bottom of the film upward. In this case, 
the depletion edge advances downward with increasing jv^j with the 
edge parallel to the film surface. It seemed plausible, therefore, 
that with Qssi^O an electron accumulation layer would exist at the 
SSI and that this accumulation layer would be gradually reduced 
across the SSI as the gate voltage, I V G J , increased. In the 2-D self- 
cohsistent depletion analysis (SCDA) discussed in Chapter 4, this 
type of behavior can be modeled { however, the program converges 
slowly when the gate field penetrates significantly into the sapphire 
substrate. Thus, it was not practical to use the program to calculate 
C-V and I-V curves and to change parameters such as the film thick- 
ness, geometry etc. 

It was concluded that a simplified program might be useful, 
especially after through-depletion occured. In this program, the 
film parameters could be changed freely, e,g. Q ssi was a P°l nt of prime 
interest. It was already known that large values of Q gs ^ caused a 
large offset in the C-V curves which had not been observed experimentally. 
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Meyer^ had observed small offsets and Heiman® had observed a lack of 
"through-depletion" (very large offsets) for deep depletion transistors. 
This behavior had been explained in terms of a layer of fixed charge 
within a high-resistance glassy interlayer watween the silicon film 
and the sapphire substrate. 9»10 The fixed charge resulted from 
ionized impurities, which might be unintentional contaminants or 
(presumably) doping impurities introduced laterally through the 
layer from the P+ and N+ regions. If doping impurities are intro- 
duced, note that acceptors from the P + region would reduce 0 ss i at the 
inner edge of the annular depletion capacitance and donors from the 
outer N+ region would increase Qssi* Therefore, the distribution for 
Qssi would generally increase from the inner radius to the outer 
radius. The simplified program allows the treatment of Q sa j. as a non- 
uniform distribution on the SSI. Such a phenomena could possibly 
explain the lack of a sharp drop in the C-V curves when through-depletion 
occurs. 

The algorithm for the simplified model is as follows? 

1. Starting with an assumed set of boundary conditions, 
typically with ip = 0 on the si-sapph interfaces solve for 
the potential field in the sapphire using a successive- 
over-relaxation procedure (SOR). 

2. Calculate the electric displacement field, D, on the 
sapphire side of the interface. Calculate the D field 
on the silicon side of the interface with the assumed 
surface potential using the one-dimensional formulas in 
Chapter 2, 

3. If D sa “ D si y % si at r = on the interface, assume 

that at this point the surface is depleted of free carriers 

17 


and use (2-17) for boundary condition at this point for a 
new SOR solution. Test each point before starting a new 
solution. 

4. Repeat (1) through (3) until no additional points on the 
interface satisfy the inequality of (3) . 

5, Increment the boundary value, Vq, and repeat. 

Additional criteria are included in the algorithm to account for 
inversion. For example, when applying (2-17), the surface potential, 
i|/ xo , tested to determine if t|) XQ < Vj - 0,55 - (kT/q) S.n(K^(l)/nj_) 
where Ng(l) is the doping of the silicon film at x a t QXt If the 
inequality is satisfied, then ij> is set equal to the right hand side 
and the B. C. given by (2-17) is expressed in terms of iJj . 

3.1 The Simplified Analysis Program 

The program listing Is given in Appendix A, It consists of a 
main program for implementing the algorithm and three subroutines. One 
subroutine establishes the lifetime and doping distributions to be used 
In the calculation. A second subroutine explicitly specifies the steps 
used in the SOR calculation and tests the residual to determine the 
maximum Value during one sweep through the field. The third subroutine 
outputs the field data for the sapphire region. The main program pro- 
ceeds with, first, a calculation of the capacitances and junction current 
vs. the gate voltage before the through-depletion occurs. The results 
from these calculations allow the establishment of the validity of 
assumed lifetime and doping distributions for the model as well as the 
value for the S1-S102 interface charge Qss2. 

After the one-dimensional calculations, the gate voltage Is in- 
cremented to a value sufficient for the gate field to penetrate the 
IS 


sapphire. Then the two-dimensional problem In the sapphire is solved 
repetitively for a sequence of gate voltages in order to determine the 
gate electrode charge and the incremental capacitance. The SOR solution 
is based on the Gauss 1 law principle applied at each grid point within 
a cylindrical region. This formulation is set forth in a general form 
in Chapter 4 where the dielectric discontinuities and variable grid 
spacing are considered. 

3.2 Calculated Results for Capacitance and Junction Currents. 

Figures 3*1 and 3,2 show the calculated results for the case where 
through-depletion is not acheived. Circles represent Kranser's data.^ 

Such errors as exist result from manipulating the data since the pro- 
gram for conditions before through-depletion is based upon a self- 
consistent 1 - D model. The calculated capacitance included parasitics 
from the stray capacitance of the probe which is incompletely guarded 
and an overlap capacitance between gate and the P+ region. The probe 
stray was measured and is approximately 1 pF. The gate - P + overlap 
at the S 1 -S 102 interface can be estimated from the intentional overlap 
(0.75 mils) designed into the structure, but the overlap at the gate 
P + edge on the SSI must be estimated. This overlap was adjusted to 
give the correct minimum capacitance. Scott etal^ have discussed this 
overlap which is troublesome for fabricating SOS integrated, circuits, 
since sufficient overlap produces a drain-to-souree short circuit. 

Figure 3,3 shows the calculated capacitance vs. gate bias for an 
entire range of Vg and two values of Vj t The dotted curve shows the 
results for a low, uniform value Q ss j/q = 7.5 X 10^ cnT?- and the 
dashed curve shows the result for a larger uniform value OsslM = 1.25X 10^* 
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The solid curve is for an exponentially distributed Q ss i with Q ss £/q = 

7.5 X 10® on the inner radius and Qssi/q - 1.25 X 10^ on the outer 
radius of the annular depletion capacitance. Obviously, a low uniform 
Qssl gives a sharp drop in G when through-depletion occurs while a 
large value of Q ss i may give a noticeable offset (threshold for Vq) 

or may even preclude through-depletion before the inversion at the 
Si-*Si02 interface occurs. The gradual reduction of C due to the gradual 
removal of the accumulation charge on the SSI due to Q ss i can be 
explained by a non-uniform distribution of Q ss i on the SSI. 

Figures 3.4 and 3.5 show the influence of the charge Q s3 2 ( a h the 
Si~SiC >2 interface) and the oxide thickness on the C— V curves. Figure 

3.6 shows the increase in the junction current caused by generation on 
the SSI as the accumulation layer is reduced. The depletion model 
augmented by a surface generation model on the SSI can account for 
increases in current up to the point where the accumulation layer on 
the SSI is fully depleted. The surface recombination velocity used 
for the results in Figure 3.6 corresponds to an average lifetime of 

2 n sec for a 0.1 micron layer. Experimentally, It is observed that 
the current continues to increase until inversion sets in, at which 
point the current falls off. The calculated field patterns discussed 
in the next section give some insight into what may be happening in 
the inversion region. 

3.3 Calculated Field Patterns 

Figure 3.7 shows equipotential surfaces in the sapphire for a 
through depleted condition while Figure 3.8 shows the normal electric 
flux density for a through-depleted and an inverted condition. In 
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Figure 3.9 the radial field profiles along the silicon-sapphire interface 
are shown for three values of gate voltage- Note that the minimum field 
location shifts toxv'ard the P + region with increasing gate voltage. Collected 
hole current flows along the Si-Si 02 interface while electron current flows 
on the SSI. At the center of the field, the electron current is approxi- 
mately 10,^, and if this current were confined in a layer 0.1 micron thick 
and the field value were 500V/Cm, the required electron density would be 
10^2 Cm“^. This value would result in a channel resistance of approxi- 
mately 100 Megohm; however, even this low a carrier density is significant 
in terms of recombination, since hole densities greater than 10 Cm J would 
result in net recombination. Apparently, all of the thermally generated 
carriers are not collected, and these tay very well be those generated 
near the SSI. After the SSI is fairly well depleted, further increase in 
the gate voltage increases the radial field and results in an increase 
in the collection of the thermally generated carriers. However, when inver- 
sion sets in, the radial field drops substantially. This must be accompanied 
by an increase in the electron concentration on the SSI. Furthermore, the 
buildup of the hole concentration on the Si-SiO^ interface will result in 
holes back diffusing toward the SSI. Therefore, it is quite plausible that 
this effect results in a decrease in the collected current with the onset 
of inversion. Figure 3.10 shows the integrated transit time from the P + to 
N + region along the SSI. Note the increase in this transit time as inver- 
sion sets in. Elementary arguments show that the total electron charge on 
the SSI must therefore increase accordingly or that the current must drop. 

The experimental results indicate a small drop in the current; therefore, 
the electron density must increase. The increase in the hole density due 
to back diffusion could then account for a reduction in the net thermal 

generation in the SSI region xuhich results in a lower current. 
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Figure 3.9 Radial Field Profiles along Silicon 
Sapphire Interface. 
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Figure 3.10 Radial Transit Time from N -N to 
P + -N Junction. 




3.4 Conclusions 


The following conclusions were drawn from the results presented in 
this chapter. 

1. The falling part of the C - V curve can be described with a 1-D 
depletion model before through-depletion and 2-D model afterward. 
Adjusting the intensity and distribution of the fixed charge 
Qssl allows the observation of three types of behavior: 

(a) a sharp drop in C at through-depletion if Q ss i is small 

(b) An offset, or threshold for |V G | , at through-depletion 
if Q ss i is large. The offset may be large enough to prevent 
removal of the electron accumulation layer before inversion at 
the Si-Si (>2 interface occurs. 

(c) A continuous reduction of C, as ovserved by Kranser, may 

be obtained from a non-uniform distribution of Q ss j. This type 
of distribution is plausibly the result of impurities diffusing 
laterally in the SSI from the P* and N + regions. 

2. The minimum capacitance is dominated by the capacitances of the 
sapphire substrate, the gate - P+ overlaps and the parasitic 
fields of the Incompletely guarded probe and lead wire. All 

of these contributions are approximately equal. 

3. The silicon film thickness has a small effect on the minimum 
capacitance observed. Small variations (±20%) in producing the 
the thin film would be difficult to detect. 

4. The effect of the oxide thickness on the minimum capacitance is 
not large but is more significant than that of the silicon film. 
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5. The concept of a surface generation current on the depleted silicon- 
sapphire interface is useful in conjunction with the 1-D depletion 
model for calculating the I-V characteristic. The depletion model 
predicts too rapid a reduction of carrier density with, increasing 
gate voltage near the SSI; hence the generation component increases 
too rapidly unless the lifetime is made artiflcally high. (See 

the level-off in Figure 2.3). Treatment of the reduction of the 
electron accumulation in terms of a surface layer for the final 
stage of depletion leads logically to the use of a surface 
generation velocity on the SSI. 

6. When inversion occurs, the radial field component along the silicon- 
sapphire interface is substantially reduced. The radial transit 
time increases substantially; consequently there is an increase 

in the electron density on the SSI and a reduction in current. The 
reduction in current is probably due to the decrease in net thermal 
generation near the SSI because of the Increase in the electron 
density and of the hole density due to holes diffusing back from 
the high concentration region on the inverted Si-Si 02 surface. 

Otherwise it was found that the modified program was efficient for 
economical generation of data and could easily be checked for convergence. 

It was found that for sapphire thicknesses in excess of 3 mils, the thick- 
ness or the boundary conditions, i.e. either zero normal field or zero 
potential, had little effect on the field in the film or on the calculated 
capacitance, current, etc. 
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CHAPTER 4 


COMPLETE TWO-DIMENSIONAL ELECTROSTATIC 
ANALYSIS OF THE SOS-GCD 

A more accurate cross-section diagram of Kranger’s 5 * 6 devices is 
shown In Figure 4.1 which shows that the gate electrode does not com- 
pletely cover the annular N-region between the P+ and N4- diffusions. 

One would anticipate fringing effects because of the lack of overlap 
on the N4* edge, and that the fringing will be more significant as the 
silicon film is depleted to the sapphire substrate. A two-dimensional 
model within the film itself is necessary to account for the fringing and 
to investigate the transit time for the movement of thermally generated 
carriers from generation to sink sites. 

Significant two-dimensional effects fall Into two categories. 

(1) The effect of electrostatic field upon the determination of the gate 
capacitance is essentially a one-dimensional phenomena until the film is 
depleted through to the sapphire substrate. After through-depletion, 
the field has a strong two-dimensional character. (2) The small two- 
dimensional effect, as noted by a small radial field component in the 
depletion region, before the through depletion will have a significant 
effect upon the reverse current-gate bias relationship. Analysis of the 
problem must be done using numerical techniques. First, the appropriate 
equations must be formulated in a finite difference scheme, and then a 
suitable grid must be chosen to span the field of interest allowing the 
necessary resolution in critical regions. 
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Figure 4,1 Cross-section of Gate-Controlled-Diode 
in SOS used by Kranzer. 6 

4.1 Electric Equations 

It was assumed that the carrier velocity depended linearly upon the 
electric field, an assumption of dubious validity near the P+ ft edge but 
reasonable otherwise. The current densities, J p for holes and J n for 
electrons are: 

J P = “ dDpVp - qp^p^ 4.1 i 

J n = qD^Vn - qp^nViJi 4,2 

The Shockley-Read recombination model was assumed so that the net 
recombination rate U is given by: 

U = (pn - n 1 2 y[r p (n + n ± ) + T n (p + n.^3 4.3 
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with t = T_ - i , a function of z. Defining normalized fluxes F^ and F n 


P n 

and generation rates G and by: 


F A Vp + (Vi|j/V T )p 


F A Vn - (Vif»/V T )n 
n i 


4.4 

4.5 


where V T - kT/q , 


G A - G/D 
P P 


4.6 


G A - 0/D„ 
n - n 


4.7 


We then write the continuity equations as: 


V • (Vp + (Vip/V T )p) « G 


V • (Vn - (V^/V T )n) =“G n 


4.8 

4.9 


These equations must be solved simultaneously with Poisson's equation: 


V • (eViJj) = - qCN^ + p - n) 


4.10 


The algorithm for a solution is based upon Gauss' theorem: 


4 ) F • nds = 
s 


V • Fdv 


4.11 


Equation 4.11 is applied to each volume element in the grid. The 
appropriate volume element is a cylindrical ring of radius r^, thickness 
Ar^ and height h. . Normal fluxes, Fg, bottom, F^, right, F^,, top, 


and F_ , left, at the elemental surfaces are expressed in finite 

Lt 


difference form for holes as: 


- F r 


P i.i - p i.i-i . hd : h.j-i tP i.i* P i,3-i ) 


4.12 


h. , 
3-1 


Vj-i 
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N 


4.13 


p. = -jui — ( p i,.i * p t-i,i ) 


p. . - p 

L Ar 


I-l 
P, , - P. 


V T flr i-l 


F = ■ + JLlL-1 i+itl _ F t.j. + P l^l.-1^ 

1 Ar 4 V T Ar 4 2 


4.14 


p = -Ui, fit.l-H , f , jij ~ ^ 1 , 1+1 <p i t .i + p i,i-fi ) 

T u 


3+1 

Applying Gauss' theorem: 


V h 
T n j+1 


4.15 


- Cr i+1 " r i-l ) (t i+l + 2r i + r i-i> . r f + 

* • — — — - f b + \ 


< r * + r^-,) Cr, 

, t 4- T1 

J R 


+ 2n -a : am h F + . “i+i z r t-i> (r m + 2r i + r i-i> 


(r 


“TT 


i+1 


~ r i-l^ ^i+l + 2r i + r i-i^ 


h.G 
.1 P 


4.16 


Defining: 


ACR(I) 

A 

( r 1+ i _ r i-i^ r i+i + 2r j + 

4.17 

BCR (I) 

A 

<r i + r i-i )/(r i - r i.P 

4.18 

DCR(I) 

A 

Cr i+i + r i )/(r i+i - V 

4.19 

DCR{I) 

A 

ACR(I) 

4.20 


and , 


AP = ACR(I)*H~ 3 *(1 - ) 

2v t 


4.21 


’I'*. A ~ ^ 1 - 

BP = BCR(I)*H* (1 -■ ’ J ) 

2Vt 


4.22 
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dp = ncR(i)*n*(i - h jl S . J j.tlrJ, ) 

2V^, 


ep = ecr<i)*h" 1 *(x - 

2Vnt 


4.23 


4.24 


where H = h and H 1 = 1/h. Similar expressions are found for 
electron concentrations except that the drift term containing ip has 
the opposite sign. These coefficients are AO, BO, Do and EO. Then: 


CP = AO + BO 4- DO + EO 4.25 

CO - AP 4- BP 4- DP 4- EP 4,26 

Finally, 

- AP*P. . , - BP*P, . , 4- CP*P„ , 
i,j-l i-l,j i,j 

- DP*P, .. . - EP*P . ... = -ACR(I)*H*G 4.27 

i"rX,j l,j"r± p 

A similar expression is found for the electron concentration n. 

In further work, a concession could be made to the nonlinearity 

of the drift velocity by testing the term AtJi/2V.j. in calculating AP, BP 

etc. If this quantity is greater than a saturation value, v *Ar/2V„p, 

sat r 

etc. it could set equal to the saturation value with v taken as 

S Q u 

8.5 x 10 6 cm/sec and 4.4 x 10® cm/sec for holes and electrons 
respectively. Values of 5.2 and 13 cm 2 /sec for the diffusion constants 
were bulk values, thus too high, and should be modified for further work. 

Boundary conditions used were zero normal current at the Si-Si0 2 
and Si-Saph interfaces and J nr = 0 and p = 0 at the left edge and 
J ~ 0 and n - Nj at the right edge. 
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For the space charge equation 4.10, we proceed in a sitnilar manner 
with additional defined terms which allow handling the discontinuity in , 
the permittivity: 

AZ(J) * 

B2(J) «= 

DZ(J) * 

EZ(J) - 

and , 

Q. . = ACR(I)*H x q<N, + p - n) + ACR(I)*0 . , 4.32 

i » J e 

where is a surface charge distribution such as Q ss j and Q ss 2 on 

the interfaces. Then, 


‘j-X /(Z j " Z 1-l' 


( Wi " Z ]-i > + “ z j-i ))/2 


BE (J) 


E j+l /CZ j+l “ V 


4.28 

4.29 

4.30 

4.31 



ACR(I)* AZ(I) etc. 


So that we have 


4.33 


“ A i,j^i,j-1 ~ V/i-l.j + 


~ D ii j ~ E i»j^i»j+1 Q iJ 


4.34 


where 


C, . - A. , + B, . + D. . + E. . 

i.j i,.1 M *»j i » j 


The boundary conditions applied are best illustrated in the 
drawing of the grid system shown in Figure 4.2, 
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4.2 The Grid System 

Choice of a reasonable grid system is the most challenging taslc in 
the analysis. One needs a fine grid within the silicon film to determine 
the small fields which exist before through depletion. However, if this 
is extended into the sapphire, an unreasonably large number of grid points 
is required to span a sufficiently large region in the sapphire. The 
accuracy of the field determination in the sapphire influences greatly 
that of the field in the film, particularly after through depletion. 

Figure 4.2 illustrates how the problem was approached. A 92 x 56 point 
+ 8 x 10 point mesh was established to cover the films and air space over 
the naked oxide and a thin region within the sapphire. Then a 10 x 40 
points mesh was selected to cover the sapphire. A linking mesh of 10 x 40 
points was chosen for interpolation between interior points of the meshes 
and their boundaries. The respective field are U + UA, V, and W. 

The radial spacing in the U field is non-uniform to allow a finer 
mesh near the P + N junction. The Z spacings in the V and W fields .are 
non-uniform with a finer mesh near the upper boundaries. An exponential 
type spacing is used. 


4.3 The .algorithm 

The system of equations is solved using the Gauss-Seidel method with 
a relaxation parameter w. At the end of the nth iteration, the field 
and carrier concentrations have been estimated at all field points. 

The (n + l)th estimate for , for example, is obtained ass 





n 


E, > + 

i.j 1,3+1 




4.35 


35 


= <1 - a)> 0* 1 + at i 

i *3 l»j i,.1 


where 0 tu < 2 and typically to = 1.7, 


The equations for P^ .. and ^ may be solved in a similar manner. 

The indicated sweep is from the left bottom corner, across and up through 


the field to the top right hand corner. After a sweep is completed, the 


order for the next sweep can be chosen to be different. The sweep order 


simply depends upon the start and end indices on two nested DO loops, and 


changing these corresponds to switching the superscripts n and n+i in 


(4.35) . 


After a sweep through the field, all of the coefficients, A. ., 

- 1 -* J 

AP. , AO etc., G, . and 0, , are calculated in terms of present 

i»j 'i»j 

estimates. Then the next sweep of the relaxation procedure is carried 


out. The iterative process is terminated when the selected residuals 


40, Ap or An are within the bounds prescribed. 


4.4 Discussion of Program Operation 


It was realized that the computational problem was very large. 


Initially the program was developed so that Q . was calculated using 

i»j 


a depletion model, i.e. ^ ~ qN^ for 0^ ^ < kT AnCN^/n^/q and 
Pl j - 0- Otherwise the potential was limited at kT fin(N ( ,/n i )/q. For 


the depletion model, approximately 8 minutes of CPU time (UNIVAC 1106) 


was required for satisfactory convergence without through -depletion. The 


procedure required 800 iterations and resulted in a potential residual on 


the order of lmV. When through depletion occurred, the residual was on 


the order of 15 mV and decreasing slowly after this number of iterations. 


This convergence was poor for calculation of the Incremental capacitance 





since the differences in two field distributions were involved. For 
obtaining semiquantitative descriptions of the field, the results were 
adequate. 

It was easy to include the continuity equations into the program. 

The relaxation procedure for p and n was carried out prior to calculating 

Q. , which utilized these values. The program listing in Appendix B 
1 > J 

includes this feature. After computation of Q., ^ , the relaxation of the 
potential tj/ was carried out. The results from this procedure were neither 
surprising nor encouraging. After 20 minutes of computation time, the 
convergence was poor. Qualitative features of the fields were reasonable; 
however, from a quantitative standpoint it was obvious that perhaps l hour 
of computational time would be needed to establish an accurate solution 
for one set of boundary conditions. Therefore, use of the program for 
modeling to calculate C-V and I-V curves for different doping and lifetime 
distributions was out of the question. 

Subsequently the program was modified leaving out the continuity 
equation. The space charge was calculated as discussed above, after an 
alternative scheme was tried and ruled out. In the alternative scheme 
the space charge was calculated using an exponential formula involving 
the potential. This scheme discussed in a preceding report. 4 This 
method increased the computing time too much for a dubious contribution 
to the accuracy. 

After the fields are calculated using the depletion approximation 


for the space charge, the charge on the gate electrode, Qgj s in the 

substrate, Q c , and the current, I CTT , are computed from sums approxi- 
blL 

mating the appropriate integrals; 
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I SIL 


D • n ds 

f 

J V 

j q G dv 

h 


A. 38 
4.39 


where S 1 is the surface of the Sl-SiO^ Interface and V is the volume 
of the n-zegion. 

Note that in using the depletion approximation G ~ nj/2r and no 
recombination in the space charge region is allowed. This causes some 
difficulties as we shall see later. 


4.5 Discussion of computed Results 

In Figure 4.3 the concentration and potential profiles along a 
cylindrical surface approximately midway between the P + N and NN+ 
junctions are presented. The convergence of the calculation was poor, 
showing inconsistencies in the radial field distribution for the concen- 
trations given by the solution. Nevertheless, the qualitative results 
are interesting. The gate voltage is -10V which corresponds experimentally 
to the condition for through depletion, although through-depletion is not 
shown here. The solution shows a significant pile-up of carriers near 
each interface, with the pile-up of electrons at the Si-sapph. interface 
being of particular interest. In this region the lifetime is very low in 
the model used, and while the hole concentration is low, the pn product is 
high enough for significant recombination to occur. The calculated 
current for this condition is approximately 14nA in fair agreement with 
the 18nA measured by Kranzer, and certainly dramatically better than 
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succeeding results using the depletion model. The location of the zero 
space charge surface is in fair agreement with the depletion edge found 
later using the depletion approximation. 


After the program was modified deleting the continuity equations, 
it was run for a series of gate voltage values and for different values 
of surface charge, Q gg ^. Figures 4.4 and 4.5 show the depletion edges 
for two different values of Q gg ^, £he si-sapph. interface. Figure 4,4 
shows that, with Q gg ^ = 0, the field from the sapphire due to the P+ N 
junction penetrates the film. Qualitatively this 3s what one expects, with 
a significant penetration into the lightly doped silicon near the inter- 
face (2.6 x 10 ll4 cm 3 doping at the interface). One should keep in mind 
that the horizontal scale is over 150 times the vertical scale when 

interpreting the graph. Figure 4.5 shows that for Q /q = 2,67 x 10 1 1 

ss 

the depletion edge remains parallel to the Si-Si0 2 interface advancing 
into the film. 

The calculated capacitances for the results in Figures 4.4 and 4,5 
did not agree with the qualitative pictures of the depletion edge given. 

The convergence of the potential calculation was suspect and the number 
of iterations was increased to 800 per point where the starting solution 
taken was the solution for the last gate bias value. In Figure 4.6 the 
left hand curves were obtained in this way. The agreement of curve B 
with Kranzer's data was encouraging. Circles indicating his data wore 
obtained from his normalized capacitance curve. 6 using the maximum 
calculated value of 13.9 pF. Curve A, with the larger value of Q gg ^, 
shows the same offset that was observed in the simplified two-dimensional 
analysis. The disagreement at the left knee was attributed to the crude 
description of the doping profile and the 0 _ value chosen to compensate 

SSfa 
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this. Convergence problems were suspected because the calculated 
curve leveled off above the data, since this was observed before. A 
determination of the capacitance for - 15 < < - 14 V then indicated 

a value of less than 3 pF with the capacitance falling with the number 
of iterations. This fell well below the data. Then the solution was 
started again at V„ = — 17V where the capacitance curve should rise 
again. The calculation was run for 1800 iterations before a change was 
made in V to determine the capacitance. Then V_ was incremented 

L? G 

in one volt steps and 800 iterations were made per point. The results 

were perplexing since the beginning capacitance value was 1.5 pF, far 

below the observed data. Otherwise the curve was shifted to the right 

of the data which could be accounted for because of improper values 

for Q ss 2 » Q ss i> N b (x), t g £, etc. Curve C in Figure 4.6 shows these 

results. At this point it had not been discovered that the capacitance 

calculation had not allowed for a significant overlap capacitance between 

the gate and p+ junction. The low value, 1.5 pF, agrees well with the 

value obtained with the simplified model without overlap. 

There was a gross discrepancy between the calculated current , shown 
in Figure 4.7, and the current measured by Kran 2 er. 5 A major part of 
the difficulty was the lifetime model which was an exponential extrapo- 
lation of Kranzer’s curve 5 which should be a reasonable description near 
the Si-Si0 2 interface. The "effective lifetime" fitting the I-V data to 
full depletion which was used in Chapter 3 had not been employed here. 
However, the results for the low 0 ^ value also show the current at low 
gate bias is larger than the saturation value with a high Q ss ^* T ^ ie fie l d 
pattern showed that the Si-sapph interface was depleted at low gate bias. 
This points out another difficultv with the depletion model which was 
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anticipated in the previous section. Namely, this model assumes all 
current in the depletion region is collected; however, even if the 
sl-sapph interface is depleted, at low gate bias the radial field 
component is so low that the generated carriers could not be collected 
before they recombine. Furthermore, most of the increased generation 
current for the smaller Q sg ^ value represents holes which would be 
generated far removed from the P+ N junction thus making collection even 
more unlikely. Otherwise, the current vs. gate bias curves show that 
the current increase after through depletion (~ 7nA) observed by Kranzer 
is too large to be accounted for by extension of the depletion region 


further into the silicon not under the gate. 

Finally, Figures 4.8 and 4.9 show the field distributions in the 
silicon after through depletion occurs. Values for the potential along 
the interface agree well enough with those obtained with the simplified 
model considering the differences in the parameters used. The results 
in Figure 4.8 are for V_ = - 18V which is before inversion sets in 
according to this model, while in Figure 4.9 V^ = - 24V and almost 
all of the Si-Si0 2 interface under the electrode is inverted. The 
results in Figure 4.9 agree qualitatively with what one expects to 
happen when Inversion sets in. The radial component of the field in the 
top layer of the silicon film is reduced and the collection efficiency 
drops. This result can of course explain the fall off in tne junction 
current after inversion begins . 

4.6 Conclusions 

The results obtained with the models discussed in this chapter 
show the following: 
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1. The actual charge distribution is probably significantly 
different from the depletion model. Consequently the 
depletion of carriers from the silicon-sapphire interface 
proceeds more slowly than predicted by the depletion model. 
There is a pile-up of holes at the Si-SiO^ interface and of 
electrons at the Sl-Sapph interface for an N-type film, 

2. The depletion model predicts the advancement of a parallel 
depletion edge from the Si-Si0 2 interface; however, uiless a 
fixed positively charged layer is assumed to exist at the 
Si-Sapph interface, a depletion occurs also from the bottom 
interface. The field strength in the bottom depletion region 
is low. Thus if such a region did In fact exist, thermally 
generated carriers in the region would not be collected. If 
one assumes that all the carriers are collected, the junction 
current calculated using an extrapolated lifetime profile is 
much in excess of the measured values. This latter difficulty 
is to some extent avoided by using a lifetime profile such 

as used in Chapter 3. 

3. The field pattern obtained after through-depletion and before 
inversion shows that the collection efficiency of the structure 
should be significantly enhanced over that obtained before 
through depletion occurs, 

4. The field pattern obtained after inversion at the Si-S10 2 
interface has begun shows that the collection efficiency should 
drop significantly with inversion. 

5. The simulations gave current vg, gate bias results grossly 
different from the measurements; however, they indicated that 


47 


the increase in the current with gate bias after through 
depletion is not simply the result of extension of the 
depletion layer into the region not under the gate 
electrode, although this phenomena is also involved (see 
Figure (4.8)). 

Experience with the program operation led to the following 
conclusions: 

6. The extended two-dimensional depletion model converges 
very slowly even though an attempt was made to minimize 
the number of points required to span the field of 
interest by using coarse, fine, and interpolating fields. 

Such a simulation model is expensive to implement for 
studying the effects of parameters on the C-V and 1-V 
curves. It is estimated that a simulation with adequate 
convergence would cost approximately $1 K at external computer 
time rates at MSU ($360/hr). Application of further pro- 
gramming and numerical analysis expertise could possibly 
reduce the cost. 

7. The two-dimensional model including the continuity equations 
is even more expensive to implement. However, this model 
may very well be more useful in gaining insight into the 
physical phenomena. A disadvantage, other than the extensive 
amount of computing time, is that mobility parameters are 
needed, and the scattering phenomena itself is complex in the 
thin film. Perhaps for the reverse bias condition the 
scattering phenomena though important is still secondary. 


The few results obtained here show that the very small 
lifetime near the Sl-Sapph interface need not necessarily 


manifest itself in as large a current 
by a depletion model. Very possibly, 
recombination may occur thus reducing 
of carriers collected. 


i. 


value as predicted 
a region of net 
the actual number 
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CHAPTER 5 


CONCLUSIONS 


Concerning the capacitance vs. gate bias curves measured by 
Kranzer,^ the following conclusions are made: 

1. The initial decrease of the capacitance can be described well 
by a one-dimensional depletion model. 

2. As the depletion extends to the silicon-sapphire interface, 
the majority carrier concentration there is expelled by 
increasing the gate bias. Three types of behavior are 
possible: 

(a) Abrupt expulsion and an abrupt drop in the C-V curve. 

(b) A gradual expulsion and gradual reduction of the C-V 
curve as observed by Kranzer. 

(c) A threshold may be observed, in some cases resulting 
in inversion at the Si-SiOg interface before depletion 
occurs . 

The behavior may be explained by the presence of a non- 
uniform fixed charge distribution on the silicon-sapphire 
interface and by the fact that the depletion occurs as a 
distributed field effect in which carrier diffusion, 
scattering, generation and recombination of carriers are 
involved. 


3. The flat part of the C-V curve exhibits a minimum capaci- 
tance which is determined mainly by the capacitance of the 
sapphire, by the overlap capacitance, and by the parasitic 
in the measurement which cannot be made zero even with 


careful guarding. The overlap and parasitic are dominant 
with 1 - 1.5 pF due to the sapphire. 

h. The rise in the C-V curve with further increase in bias is due 
to inversion of the Si-SiO^ interface region. 

Concerning the junction current vs. gate bias curves measured by 
Kranzer,^ the following conclusions are made: 

5. The beginning part of the curve rises because of collection 
of carriers generated in a region which is continually 
expanded as the gate bias increases. This portion of the 
I-V curve is useful for obtaining an effective lifetime 
using a one-dimensional analysis. 

6. Depletion through to the silicon-sapphire interface precedes 
more slowly than predicted by a one-dimensional model. 

However, it appears reasonable to treat the depletion as 

a surface phenomena with an increase in collected current 
due to a surface generation current as the depleted area 
grows. After depletion of the interface there is a small 
(12$) increase in current which may be attributed in part 
to the decrease of the radial transit time and the lateral 
increase ( 6 %) of the depleted volume as the gate bias is 
increased further. 

7. When inversion of the Si-SiO^ interface occurs there is a 
substantial decrease in the radial field and an increase in 
the radial transit time, which should be accompanied by an 
increase in the electron concentration on the silicon- 
sapphire interface and in the hole concentration on the 
Si-SiO^ interface. These increases in carrier densities 
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should result in a reduction in net generation; hence, the 
collected current drops as observed experimentally. 

Concerning the utility of the Gate-Controlled-Diode for lifetime measure- 
ments, the following conclusions are made: 

8. An "effective-lifetime" profile can be constructed from the 
C-V and I-V curves using the depletion model. This profile 
will be meaningful for data obtained for, probably, 50 - 75 % 
of through-depletion and can be extrapolated to through- 
depleti on. 

9. The actual generation lifetime will be less than the "effective- 
lifetime." By obtaining the current at through-depletion and at 
saturation for a large Junction bias voltage, the ratio of the 
former to the latter is obtained. The "effective lifetime" may 
by multiplied by this ratio to obtain a better estimate of the 
generation lifetime. 

Concerning the utility of the simulation programs, the following 
conclusions are made: 

10. The two-dimensional programs are too expensive to use for routine 
interpretation of data. The simplified two-dimensional program 

is useful for estimating the minimum capacitance and is economical. 

11. The two-dimensional program including the continuity eouations 
appears to be the most useful for obtaining insight into charge 
distributions and carrier flow, but it is very expensive to 
operate. Further refinements and operation of this program 
should be carried out. 

12. A more extensive one-dimensional model may very well be useful in 

studying the falling part of the C-V curve. It would not be 
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predictive of current flow but would probably require input data 
from the current measurements. Nevertheless, more accurate 
modeling of charge distributions may be obtained using a model 
including the continuity equation. Such a model vhould be 
inexpensive to implement. 
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TWO DIMENSIONAL ANALYSIS INCLUDING 
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■ irtEuSION a* 25»U) #HA< 13,7)*.^ f 93*56) 

iUMEhSION AAR126) »6A|U26> *UAR 126) * EAR(26> 
niMCiviSiOM AAZ<12) »UAZ(12) *PA?ii?) *EA7< l?) 

r i. t'-'i F hS IO M ACR ( 94 ) # bCh ( 94 ) * L»CR 5 9" > ♦ ECH 1 94 > 
niivCuSiON hCZU> 1 ) ibCZ(bl) *DCZ(6) ) *ECZ(6l > 
niWE^SlON a< 100) 'CS(bl) »°HI (51 >»TAUf 51) tZBUft) 
iUpCuSION tM93«51)f n(93t*n iG(9**»5I) 

PuAU Nl»NO» f 'l2 

nr=S*2 

Pu=l ./LrP 

r n=lo« 

ru=l./UM 

fu = l %5L1 0 

; ±2=r* 1**2 
i J«KI i — u 

-j'jl-O. 

-oS2=2.67Lit 
l. i OP=bu 
^- 1 *cE“ 1 9 

~:»S2=ab52*w 
Vl=0„025q 
T l-L- u • b 
T uX=u*4 
Tjiru.b 
r j'j-iy • L 1 5 

TmUO=2.F-c> 

'>uTA 1 = 0.92 
^i_I A £ = to* 3^ 
p±=3. l4lo 

fpSA=U.3+o.A54E-14 

FrSB=U.7*o*i)54E-i4 

F^bC=3.7*o.^S4E-14 

Fr SO-a.ASHL-I 11 

!*=9 

Ti_=L 7 

To=2b 

I.a=a 

L«»=a 

Kx = T u 
J*=ti 
’* t**u * 

’■^=lb2» 

W«J— 14-0* 

L 4=8 

I*t=9 

Lb=3^ 

T J-4i 
K*-=lu 
Jt-- 1 i 

?u=3b 
L/=92 
T 

10 = 8 * 
v o=bU 
Jo=5i 
K-*=5b 

Kb=6U 
,10=61 

H=TS1/4J3-^) 

Pi i=l • /H 
mj1=100«*m 
Hj2=iOU**m 


ORIGINAL PAGI3 l-j 
or POOR QLALIXg 




1 


t V ,b/VT 
T 1 # F — l> 

t mu i- 1 » r— o 

r u .*in u rL cut.FFiur.Nr-. i<, -.-field 

1 1=1 » 1 x 

; vU)-(I-li*wi/Ll 
4 !=Jtr±P 

1 1 Mi)-Wl+tl-II)*W2/L • 

f 'w r, l=T 2 n? 

'' * II ) -Wl + W«-+ < 1~12 ) *W VL3 

Ti--- L J"1 
'•u ■» i-?,l0 

.T«K( i )=PI*A( I>*{ X( 1 + D-xC I-i)) 
•^KUUPi^vm+vd.i) )/(x(d-a(i-1)> 

•/m( ij^PT+txi r+i)+xd> )/cx(i+n*xtn i 

y - n,< (I)=Pl»A(I)*tx(Ifl)-xU>i 

cl* 

in yY I )=-1SI*(4lPHA + *(.11-J>-1 # )/(AI..PHA-i* ) 

” W 1 tu j-P * M 

*.mUU-FPW(X< J)-X(j-l) ) 

' MUu)~FPbfc*(Xl J+ll-xU) 1/2. 

'-W(U-PAU J> 

IP '-mZ(o)=FPpA/{X< J+I)-X<u) 1 

"mLCuLATT LPtrucif-IS FOR V-FIFLO. 

U 1h 1-1 » 1 4 

14 ' v 1 > - ( 1-1 ) * v i/ 14 . *LA ) +3.+W1/4 • 

’ v io 1-14*15 

16 'i x )-( I-IU*wP/LR + wi 
T v= 1 b-l 

to 1 — 2 * I* 1 

■■~K(i!-PI*A<I)+ X(IU)-Xd-l) * 

Puk D=Pl*iy C T )+y 1 1 - 1 ) )/(x( I)-X(I-n ) 

• uKU)iPiny(i + i)r<( D >/(xt n-xti-D > 

1* r «-.K( i )=AbMT) 

' >_HHh — 1 « ?J 

" n ',!•(- { hLI'HA *+K2“l . — K A * C ALP M A— 1 ■ ) )*<ALPHA-t , ) 
^=(r.2-1 ) *m! PIIA*+KP-i. 2*ALPHA**tK2-l )+J , 

'« = -#Fl/AP 

m( piwt + jA 

’m r I u-kPS i /ALPHA ) 

T r ( P>t T i U U #f)*01) O' Tu *»U 

" u =TjI/ ( Jo-U 
■'w Pc. 

r? J i HO 4 l alpha** ( ,J c -j) -1 * ) / t Al.PHA-l . ) 

w i? *t U— p » h, n 

< w(u)=rp JW /(7iJt j)-/ u u-i ) ) 
oL U)-Fp o „H70( J+1)-2U{J-1 U/2. 

VA f a > — 1 <P ) 

P4 -u/(o)-Fpb«/(,7HCj+l)-2u( J}> 

'■wLFf IcTFujF FOP tJ-UELll. 

'■u Po i = l*lf> 

:f> /il)-(l-l J* , * , 2/t2.*Lr*LU + .ni 

"u Pu l=lu»T7 

PH MD-X1I6J M 1-16 )+*,?/ (c;.*L5> 

Tw=I/-l 
*w 5u 1“2 f 1' 1 

( 1 )=PI*A( I )*U< I + l)-<{ I-ll ) 

nuvU)=piMV(i i+xn-D >/(x( n-x( i-in 
-wk ( i)-PI+tX(I+l )+Xf I ) >/{X( I+1)-X(I) J 
”0 e*Vrt(i) = ACHtT) 

Jc-P) 
y j^-i j 

v 1 3 ) -n . 

-o it: ^=4*0^ 

;;2 V ^ Jj -TbT* i J-3)/ c J3“^) 

<u 5 h J=uo*J4 

74 v t J ) -X t J3 H ToX* ( J-J‘< ) / { J4-J3 > 

36 11 IJ) -X t J4l +TLL+(J-J'. J / t J r )-J4> 

lw^<^J=rH^ M /Txt2)-Xf 1) ) 

■v-2(t:>-FpbM*(X<3)-X( 1) )/->. 

Mw/ (t.)=PC^U) 

-u/( c .)=FPb«/(Xl3)-X{^) ) 

-»U< j1-FH^A/iXt3)-X(2> > 

f 'u2( J> = <FFbB*<X<U-’'(3) >+LPSA^ (Y(3)-X(p) ) )/P. 
-^Kil-FPOu/ tx(4)-x< i) ) 




61 




7 fl 


uO 


42 

43 


7r a 

a no 
ii(.2 


4 4 
'I 5 
4 6 

4 A 

50 
1 2 

1 n 


3b u=4#t\7 

Av2(J)=rpbb/(XtJ}-X<J-l) ) 
•<bZ(^)=FpJb*tXtJ+l)-A( J-U )/2. 
r 'v^(u)=PC <Ll J) 

p*2<o>=ppbo/<x<Ju>- A <o) ) 
Aw2(o3>=EPSP/<X{J3)-. A (J3-l) > 

nuf(w3!-qt2u*! X<J ' 5fi, “ X( Ji) )+t - nsP *^ x(J ^>' 

rbZ(i3)=EHbC/U<J3+l )~x< J3> ) 

Jo=J4~l 

Ju=Jo+l 

Po THi_ oXIUl. RFulON# 

~'o 4b J=JL#v'U 

.^ 2 {^)=rpbC/txtj)-xf j-i) » 

' ? u2(^)=FPbC+ {VU+n-X( J> )/2. 

n u2 ( u ) ( J ) 

^o2<«>=EPSL/(X<J+l)-X( J) ) 

Aw2<o4J=EF5C/<A( J4}- A tj4-1> ) 

J o4 ) = ( bPSO* x { J4 + 1 ) -X C J4 ) ) +tr SC ♦ ( X < J« ) • 
+C. *— * J4 / 

J4)=EH5 n /<Xl J4+t )-X( J4) ) 

Ju=J;>-l 

4< 0=JL«JU 

'o2<o)=FpSO/tXl J)-X(4“i ) ) 

'^2(u)=FPbO* (X( J+iJ-X( J) >/2. 

^<-2 (u)-PC2( J) 

Fb2(J)=FPbO/(XtJ+U- A tj) ) 
r o 43 b=3 »J3 

rot J)=CSS*tXp(-f)FTAl *{Tsi~At J > ) ) 
TM VH!r^uo*FXP( :: 0ETA2*( TSI-X< J > ) ) 
PnI(o>=VT*bFTAl+X( J) 

~ l u 7 u a Jrr3f J3 
7urt 

0 i J, f u ) *-CS ^ J J 

■n i »o)=pi**?/cs( j) 

Cof;T if JUF 

VJ* VO* W f, lT 
FoK* a MT(3FA0*3) 

Tr(Vo.bT.U) STOP 
"i.KI l=b 
\/o=Vo~i • 1 

IF < ;ti \ II . G I . \ . 0 ) MWK r 7=1 
Sr * A U L Fi AL n noUfjOAr y VAI ULS. 

44 1 = 1 n’ 

: i I * ol ) =Vo 
n« 4 b i = 1 » i 4 
11 i I * o2 > =Vo 
’ 4o i = I<£ » T .3 
■ i I »ot)=Q. 
r,| u 4o o=3 
MU »o)=Vj 
i I7» J)=ptii ( j) 

'o So o=J^»v»4 

1 : v i , u> =VJ+ U-J 3 > *( V''~ 7 j) / ( J 4 -J. 5 ) 
mp, J)=Piii(J.3) 
o ">C. i = l»I° 

’ I i 1 f o4)=Vo 
ru=i + I7»Fo' 

J " w 54 J-I»u 

’ bi ( 1 1 J ) =Vb 
m( I bf J)=PHI l J3> 


' 'j , Su 1 = 1# IP 

l»6i=Vo-{ I-i)*( V ',-PHI <J3> )/(T7- 
s l Art | IN ,,-FlFLU Win; RELAXATION* 
Ir'ASb— U 
f **RiA A =U 
IP Ob '‘vNTifJUF 

rrAS J =lPA'->S+l 
Im=U + 1 
tu=Ii.-l 
'^67 1 = I A # 1 ' l 
Tr=4-» I**3 + il 
f? "iI.oU=VUP#02> 

* 7 

^iJ=U • 
l_uOP=l 
iuu=wt-2 


Ia) 


X(J3-n ) )/ 2 . 


XU4-1 )))/2. 


ORIGINAL PAGE IS 
OP POOR QUALM 


62 


1 1 v>“ 1 

JuO — 1 

rt ( Ji_O.LTit) JLOr? 

J i_ !> - o • 

'u ho u-jLU»vJ fl i 
•"o ho L-2tl'H> 

■.v/Ld-W £ T »wjJ 

, '.-A ( i )«Am?( J) 

( 1 ) tuii 7 ( J) 

' 1 — L* A u £ 1 ) *l>h 7 £ J J 
r-u \i\ £ l ) 1 J ) 

r -A*o f i-J+F 

• IL-(A* rt li» J-D+H* (I-l»JHX)*i#<IM#J)*r*WU»J+1 ) J/F 
U »o)-(l.-0,'.|FoA)*W ^D+O'lEGA + WTTL 
I » J>-v.OLC 
• ■ bS- kCS+u h * * 2 
CwnTiNUF 
■■ £i»o)-W{ti .0 

nntJ)=rtiiHp»o) 
ouUT INUF 


“b 7u i”i»lT 


'■wNTaMuf 
I y>OP-LuF)P + i 

Tr t LuOP »Ll . 1 (J ) GO In bu 
Tr £M„KiT.fc.o.U) Go To 77 
bGKT £klG> 
olTt_(b»buU) 

'-wK''mT ♦DATA J- OK FIEUnn*#/) 

r\ITi_(b»b4£i> FMS 
r wK^'rtT t / » dOX t 1 1 » pb, ^ f / ) 

(\ITt_ ( b » 504) 

c^WT(/KUVf'AlI ( J) In ADJAC r MT FIELDS*/) 

" l '-i 7d i' f i— i * J] 

ruuTiWUF 
l\ I fu£ui bu /) 

c*w K inT£//»^nx» f NEXT KlLiO* »//> 

'y 74 P=iiJl 
J-Ot+l-M 

•’i\lTc. t o»buo) ( w£ I • J) « i~j 1 1 1?) 

C^IjTiNUF 
•;oi T l ( D * 5U 7 ) 

7o N“lfu1 

J-Jl+l-W 

■\iT[_(o»buo) ( w £ I » J) , i = i2» 1 3 > 

C^nT iUuF 
\ X Tc. £ o * bLi / ) 

FuK'-'rt T ( tQAr^FlU.3) 

OwrJTiNUF 
Hu Tt-iE V-Fi^L 0 * 

LwCAlE poUiMDARlf VALUES n rt AUL- U FIELDS AMD 
I »i — 1 + 3+L1/4 


twt If J) ,i=i2,l3) 


rt anl- u fields 


Tr;TF«POLATr # 


T a- T d-i 
Hu 7o I = I»irTK 
V't.0?-* UfKl) 

Vnl = * (I + l’M ) 
i-'-4* i I-Iw) +1 
< iK # x ) =VLO 

7 £K+x# i )=VLOP ( VHI-V1 0)/4. 

■/ v K +<£. * 1 JSVLO+ ( VHI—Vl u>/2. 

I ^K + 7>» 1 )=^LO+3.* ( VMf-VL0)/4* 

' iK + 4 * 1 > = Vhl 

^uhli.iUF 

’ ° d O— ? t %j2 

c-(7u(J)-/ba))/ZB(l) 

■/Ufw>=W(ii#Kl )-S*(p (Ii, fKl) ) 

Vk lbf ))=w£i^fKl)-S*(tt(l2»Jl>-rtlT2*Kn 

CwH’l iltUF 

Ir ( T t AbS.oi ,1) GO To Bb 
M14 + 1* J2J=H<1V»3) 

!^0=i4+? 

”v 34 

+ t T-iLO') 

' l i » ) -I I £i*i f 3 ) 

rvuTillUF 
'.buTiMUF 

■’i_LAa «tfL V-*FlLLU» 




T2#Kin 




oO 


5i.M 


bio 


9? 


c 

c 


C< 4 
51? 
'>5 


c - 6 




mo 


lr.2 


1(4 


GO 

00 


Ot> 

'5b 


IU * AtXMC r MT F'lFLOS*/) 


(v<Xf j)»isj fl si#a) 


ViLOrt = l • 7 
LuuP-1 
'uF=i5-1 
Ju.0=d2~? 

JuO-JLO-1 
Ir < JLO.t T«*0 JLU=P. 

"uS-u* 

"u ^Li Jui>t JHi 

fio l=2*i'JP 
1 ) *Mb'( JJ 
•-uRk ( J. J ( J) 

"-up t\ < i > *uL-? ( J) 

r-A+u^D+E 
VuL‘>-VtI» J) 

f i IL-U+vU* J-u m*v( 1-j » J)+D*v U+l*u>-»r*v( T r j+i > ) /r 

7* 1 »y) = ( 1 • uA) ♦ t V^L r 54rjM£( J AtVTIL 
n*z\r 1 1 ( J)-v"LP 
“i-b=i\cb+0V*^2 
Cui-jTinur 
ru.bnNu^ 

! uUP-LuOp-*-^ 

ir (Luop.Lt.ma) 

*t ( LPUT.u^.ij} 

( v i Ti_ ( u t boo ) 

^wK h>n/,2u Y f *UATA ro« vM'J) FTFLH*/) 

-'I'tS-bOuT (uLS) 

• Tu ( O » btJtt ) I' ns 

KlTL(6»‘ji.U) 

FuK VT(/»£^PV( i# j> 

"*u 1c. t*,— 1 T sj? 

J-J2-M+1 
■ *\i Tu ( O r Si*, ) 

CuiOTiUUF 
vt\ITt_ ( b * bo 7 ) 
r ‘u Oh h=l»u'’ 

AiviTu (b r bA<-> 

1 *wi«TlHb r 

cuh;v,A*r( tOA# 1 lFb.2) 

Fu^Ti'JuF 

">u Trifc: PEuA v AT 10^ R-K Th r 

LuAD VALUED F|<o f * 

iiil» 

"i Ilf U-\J i A bf J2“?> 

! ! i!7#a)=VU e i# J2»u 

ui7»l J=V(iU+l» J2-2) 

'u9A i — Iu»T7»ii 
•«-i‘» + 2+( I-aM/2 
'Ul#i>sV(m# Jl-2) 

Cu.iTAUuf 
■'u 9o 

'HiaJ:im*nPM)*(u(i?»i)-untin/ib 

/ ~uHTli iUF 
“.u tud 

nU#i)=t»U7#i) + (I-l7>^(u(33#l )-H(17»l) )/m 
OuHTiNuF 

r up : i 7- i D^Up'-TA v 

'V lud I=j*i IUP*2 b.ca-Ui,*/ 

• *a»i) = ru(i-i»i)HJti+i»D )/?. 

CvHTibbF 
Oi_S-u • u 
t. uOP-1 
Oi'.LOa^I • 7 
CuHTiNUF 

+ ♦♦ * J 

IuHma-ICmaa+1 
C«»iaX=ICMaa*CS5 


U-FItlP 

THE \/-f IE( 0 AFt> IHTF^rOLATE. 


It 6 


r >i iSFu * 

CALLHlmP THE LATES1 OEHERAT Ton RATE. 

. ]-*3 

f> w W l ob 

yl’o>-1 .b-T+ACRd) *nbl*n-jI2-P(T,J>*o(I» J) J 
*/vlAu(o)*(p( 1 1 J) +0 ( T » J J ♦P.+HX ) ) 

f’.u iu7 J=a#K3 

'w u / t=^#ihp 


T 


mu* 4 • ■ <*-« w*o i 


Tiir n m m 




64 


lr 7 
1 »:H 

L 


n7(l 

fc>7? 

674 

7nn 




*L-> J *ACP ( I ) *h*(f.iI e ;-Ptl»J) *0( I * J) )/ 

* ( \ AUt J) * (PI 1 1 Jl+0< X f J)+?.*UI ) J 
J-J3 

no t=2#tup 

G i 1 t j)-1 t b-A*ACR<I)+iiSL*UU2“PUtJU0(ItJ> )/ 

♦ ( i A'U J)*(P< T t JHo< I» J)+p>. *WI ) ) 

Cut iT^t ior 

TnK m £a] *5 Lf>jPS CALCULATE T»*F CnRRIFp COUC rt iTRATlUrj«; . 

■'ViLCjM—U • S 
J-j _ 

T up- 1 7-1 
f 'u 7u J I— TOP 
c.-ACiX < 1 ) *KH 
n-BCt> ( I ) *rl 
7-UCn ( i ) *tl 
r-t.C«(i)*Kh 

AP=*u(I.~i v*(Uii» j)-u* it j-ii ) 

Av=A*U«+Tk*(U(I» j)-u( It J-ii ) 

^= r u(l.-| kUUUrJ)-u(HrJ)) 
n u— ,l ,* ( 1 • + I V*(U( 1 1 J)-U< I-lt J) ) 
w+(U( It Ji-ud+lt J) ) 

^w-'}* ( i.+iv* (-m 1 1 J)*ui x+i» J) ) 

^r=C* ti.-i w*(Ut 1 1 J)-IJ< 1 , J+l) ) 

"v/=E+ ( 1 . + I V * L ' U » J ) -L (I t , 1 + 1 ) ) 

Cr-=Au+i-iO+JO+EO 
Ct«=At- +bP+uP+£P 
CP=O'-a0 
Cu=Ch“AP 
PwLJ-PITtoJ 
= 1 r Ji 

lr (I.Eu.luP) GO TO tn70 

•>i It-(bP*p( T-l r J)+DP*P( J+it.J)+LP*P(Tt J+1 )+PO*G{ JtJ) ) /Cp 
Tr l t.eo.l) Go TO 67 > 

1 1 IL— ( bP*0 1 T -1 t J) +Up* 0( J+l t J) +£0*0 ( T t J+I ) +fjn*G f T # J) ) /C(i 
Go To o 74 
OL = Cr-L/0 

’’ l !L-(uP*P t T -l» J) +E n *P t 1 1 J+l)+pn*G( T * J) )/CP 
'■> v i +i t J) ruo*pt 1 * J)/ V 

•>1 iL-tu^+ut i-i» J} + 0 ^+ 0 ( i + it jhlp*o( t, j+i )+tn*^( t t J ) )/rti 

G^ T J 074 

C«.-Ch-LiP 

0 i lLb(U0*otT + 1 # JJ+tf'+Ot Tt J+l )+H n *G< T t J) )/C ft 
v i-i t j)=up+0 t i » J ) /' : \j 
CotJTiNOF 

''Uiw) = tl. -o;^roA ) *P^LO+o M EOA*PT T L 
' d»o) = U«-0,«lEGA)*Oi.o+OMEGA*OTIl. 

Tr (^11* J) .LT.P* J PtjtJJcn. 

Tr (Oi It J) tuT.n. ) 0< T tJ)=0. 

TrlO(ItJ)* 0 T * CMAX ) MItj)=CMAX 

"(It J)=CMAX 


-rClL 


Lv 


Ii l^>i 1 1 J) .oT.CMAX) 

Cu.MTjJIUF 
Go 7o4 J— H#K 3 
no 7 u 2 I=2tItJP 
^ -ACts ( 1 ) *t\H 
M~lSCk ( 1 ) *h 

' 1 -bCr\U)*n . . 

F-tCn ( i ) *i\h 

Af'=A*(l»-»V+(U(I t J)-UClr J**l) ) ) 

-\u=A*Q. + l'tf + (UlI t J)-U< It J-l> ) 1 
.ir=n*U.-lV*(U(ItJ>-uU -1 ij) ) ) 

>o=T»(i. + I v + (tKIt J)-U(l -1 t J) ) ) 

u.-ntiui it j)-u( i+i t <j > > > 
w=G 4 d.f i v*di(lt J)-b< i+ Lt j) ) ) 
n -n+fl.-l w* (UC It J)-u(I, J+l) } ) 

-u=r + u. + rv*(u(itj)-uUt J+l ) ) ) 

Cr- Ao+L)0+b0+E0 
Ch=AP+uP+LP+EP 
"oLO-P ( T i o ) 
r'L.L=o( i » JJ 
lr ( I.Eu.IOP) GO TO 

Pt IL=(aP*P(I t J“1 )+Ur*P( I-lt J)+UP*P< T + lt J)+Tip*P(Tt J+l ) +prt #(3 ( T r J ) )/ 
*Ct' 

T» ( I .Eo.^J GO fO eES 

n * I L-( AG+O ( I , J-l ) +tjn*rO ( j -1 1 J) +1:0*0 < T+-1 1 J) +E0*0 ( 1 1 J+l ) +Nn*G t T , J) ) / 

♦ Gn 

Go To 0 A 4 
CP-Cr'-oO 

r ’}iL-(AP*P( T ,J-l)+Bf*P(i-l,Ji+LP*P(T»J+1 )+GU#Jl)/CP 65 

I+i. to)=uo*P( I * j)/r»> ■ 






6 i=P. 


bt’A 


70? 

7r>4 


690 


b«? 


6 ° 4 


7(i6 


C 

O 

L 


C I IL-(AO*oiTf J-l ) +h'’*o( l-l r J)+l/P*P( T + l i J) MlO^nt T* J+1 >+;^ + G< T, J) )/ 
^o To of*4 

C it-Ch'li*’ 

!L-( A0*01T» J-l ) +U^*0( J fir J)+tP*0< T » J-M ) ♦NP*G( T# J) > /r , 

‘■'il-Xf J>=ljP*0(1 » J>/> 0 

Cout xikjf 

*-Of./|F(»A) *PnLi)+0 M EOA*HTTL 
Tr(PtI»J>«LT.O.> Plr,J)rn, 

^ilf j)=(l.-PMFGA)+0Lo+0peGA*0TU, 

(0(l*J).LT.fl.) 0(TfJ)=n. 


i r 


o( I# jl-C'^AX 
Pt I t J)=CMAX 


If ( 0 ( I f J ) .OT.CMAX) 

It (pilf J) .oT.CMAX) 

Cont xmuf 

CoWTlNUF 

J—J3 

M 7u6 T=<if TUP 
A-AC.tt I>*KH 
R-UCk 1 1 ) *H 
''—OCR J I >*M 
^-EChl 1)*KM 

•;r=A*( i.-iv*{ini»J}-uU».Hi> ) i 


Ao=A*<X*+i v* (Utl r J)-uJ I i.t-l) ) j 


r t r =R*(i.-!V+(U(If J) 

->o=M*<l.+Tv*U«t I » J)-u* 1-1 • J) ) J 

n -=n*(i.-* v *tui I» J)-u<l4l» J) > ) 

•^w=n+ ( n+Tv * cm I * J) -u< i4 1» j) ) ) 

\t =E*U.-Tw*(UUf t+l) ) > 

^o=E»U. + l\/*(UlI»J)-uCl t J+I) > ) 

'V-Ao+X>0+uo+£0 

C,.-AH+bP+op+EP 

ci-cr-lo 

Cut— Ci'i“EP 
• J oLP— P ( 1 t J J 
0t-U=O( 1 r j) 

TriT.EO.lUH> GO TO 6‘iO 

PUL-i Ar*P( ft J“i) +BP*P( I — 1» J>+UP*Pt T + l » J)+PP*G< T# J) >/CP 
Jr ( I«EG*2> GO TO 69? 

r i IU=i A0*0(T» J-l >+bn*0( i-lfJ)+U0*O<T+lfJ)+ND*G( If J) t /<"r f 

Go To o°4 

rr=CH-00 

p i 1 L= t MP*P ( T » J-l ) +UP*P ( I -1 » J ) +H^*G ( T * J) ) /Cp 
Pi I+if j)roo*p( I» J)/np 

0 1 lL-< a 0*0( T f J-l ) +60*0 ( 1-1 fJ)+L)0*0( T + l fJ)+f ‘0*G( T fJ) ) /CfJ 
Go Tu b 04 
Cn-Cu— 6P 


M IL=UO*o(I# J-l ' +Uo*0 ( I + 1 » J) +UP+G ( T t J) > /CM 
o 4 I-a» J)=OH*0( '/'-j 


CohTiNUF 

Pi 1 #o> = ( l.-Of.irGA) *P^LU+0 ‘ j 1EGA*HT TL 
C 4 i * o) = ( 1 .-OMFoA) *0! U+CW6A+0T It. 
ir (PI 1 1 J) *LT.O. ) P(J»J)=n. 
ir tot If J> .LT.O. ) 0 (l»J)ro. 

' ‘ OT.CMAX) 'MI # J)— CWAX 

•UT.CMAX) P < I » J ) -CM AX 


original page ip 

op POOR QUALITY 


********** 

POINTg. 


1(0 


no 


111 


66 


tr (01 If J) 

Tr(P(IfJ> 

I'vyliTilluF 
"viCGh — 1 • 7 

EhO of O^llCLf 'TR mTIOM t OOP 

***i+******* J 

compute Tiir spAcr charge at the fjelp 

lu 9 Irtfl'jp 

Pi iOPvj+ (CG (o ) +P( I » J) — o ( I » J) ) 

^o( If U> = l.t_-ii*ACR( 1 ) *H*rHO +l*t-4*ACR< I )*0$S1 
CouTiMUF 

^o no 

^o no i=tif nip 
"rtO=o* (C5(J)+PlI»J)-ullf J) ) 

Ooilf J)=l.t-A»ACtUl>*H»RMO 
CoUT ii'lUF 
J-J3 

r o 1*1 IZ^lflUf' 

PMO=u*iCSio>+P(Tf J)-o( It .1) ) 

)j(It J)=l.t-B+rtCK<I)*H#p‘«0+nE-4*ACf5 ( I)*0S5? 

rowf *i jU^ 

'uLAx (HE U(l fJ> Flr-bD THROUGH GAPP, #SIL. f AMP ovIDf. 
JitI-o4-l 
Oo Iio Js^iJHl 
9u ll4 T^JfluP 


T 


114 

116 


110 


120 

122 


513 

514 

515 


C 

C 


124 

126 


516 


A-aCu ( 1 1 ♦hU? ( J> 

R-bCidl ) ♦uL''( J) 
p '-UCi\(i>*ut?( J> 
f--LC,\( i >*LC7(J) 
r4A+u+U+F 
t‘oL')=U( T , j) 

• i i f V3 *=?V >5So l i^i-'tp , iA*SV{J +i * J) ****•'« 1 • J* i > « t , iu/r 

U (•!( 1 t J> • vjT.PHI ( J) ) u( T » J ) — PH T ( J) 

u=n(l, j)-u^un 

CuhrAilU 1 " 

CoijT itJU^ 

Tt- h<LS,r?T.i.tl(J) RtS^l *t 10 
'm5=5'JF.T(klS/ f 1 1 7-1 ) * t J 4 -I) ) ) 

Tr {'.\,tUT.of.0> WRtr r 't 6 »» 5 U) 
ftvH 1 1— PkS 

■ U 

ri_o=iH+l 

■Iw 11 H rz:i(_ 0 ,IUP 
'J-ACk* 1 >*hC?{ J4) 

R-bC»\( 1 )*uL 7 ( J4) 
i-L>CinU)*pC7(J4) 

"-uCk < I ) *lg7< J4) 
t-a+u+U+F 
‘ VLiJ-U ( T t OH ) 

'•-l + l-ifl 

'ft lL=(A*UliP J4-J )+b+u{ i-l » J4)*0+LM 1 + 1 * J4 >+E*UA <"» 2 ) ) /C 
’ 1 1 1 1 04 ) -Ti . “O^tGA ) *s l 0LD+'iMt.GA*UTlL 
'<m( olJ=UU * J4> 

'u='Jl I » J4J -'I0LU 
'^5=r(Eb+DU*+2 
CuNTlNUF 

^LAa the UA-FIFLD. 

. it-0“u4+ 1 
Jiil=u5-1 

">v 1^2 J=JL0»JHI 
I.-J-J4+1 

'"o IcO TziLOrlUP 
"-H-1-I8 
A-ACrt { £ ) * mC7 ( Ji 

'*~aCt\ £11 *DG7 ( J1 

--EC.\U >*tL7CJJ 

C-A+u+u+E 

MwL J=UA ( M »*_ ) 

’!' ff : - t ( ^iyv iM A!7Z’ril ti * rU ^ (M " 1,L,+LjtuA(M ' f l*L) + E r *nA(»'#L4U >/r 
i!rt( 1»L)-(l,-0 /EGA) ♦UoUj+OMEGA*UTIL 
"'•u-l Jm ( M» L ) -POLL/ 

J cS=KFb4Du**2 

CuNTIMUF *e 

CuHTliJUF 

K^fTAi«5}?i52JilTi ifftan > 

K&fiilti 1 .* !S.S? ,TF ‘ 6 ** 13 ' rhs °»^ a 

'-ofr'!*T<?cUf < “in.3) 0» 

i uUP-LOOp+i 

I> (LuOF.Lt.LTOF) GO TO 184 

C '" , ?^uTt ,t T l C .F A cgl R ^? he0 IN ™ r S1 ‘ 1CWI LAr ER. 

|uP=u3-l 
Cuk— u * 

O~>IL-0 • 

^o leb J=4fJljP 

1^4 Tzdr TUP 

^jl'L=OSlLf i*E-4*QS( T t J) 

*'GK— uUK+G £ j. * <J) 
r'uf^T iHUF 
rubTiNUF 
r uK=i. .c~4*0i*CHR 
ivi T u ( o# bit > ) V J t \J($t ^’S IL» OUH 

rv/K 0,T(/r AOV, *\/J=* t r ia.3 f >, x # M/G“* »P»0*3»b¥* »GSTt.= * ♦ Pi ft, ^5** 

* f *£.1 u • 3 J 

Tr (r u Ru.U .°. ri U4) Ty *U7 
'■’U4»42>-G( 17# 3) 

T i-u^i4+2 

1j 4 r^iLbflb 07 


i :«4 

517 

536 

54 0 

537 

54 2 
5.38 

543 
5u 4 

518 

A?8 


132 

52 G 


522 

136 


Vil*u2)=Uli*i#3> 

Gu T v> 18 {)u 
r U) JT if UJ^ 

•JhITlIo* 53 b ) 

r -K'MT < / » 30* * * AHbRLVI ATE^ U-F XU 0 ' / ) 

J ’w 54 U M- i • J 4 * 5 

j- 04 -is+l 

>i\lTfc|o»l>4 *> (U(i,J),Ui f 17*0) 

<\1 TL t ti# 5 J / ) 

FoWImT* /# 50X» *Afn«RC^iAft.l> li-FitLD'/) 
nu 5 h 2 M = j # j.- 5 #b 

J-J 3 -M +1 

nlTt. tb#i>*»4 > <O<I#J)*I=i*I7*0> 

ruK'Vn/* JOY* ♦AUaREi'lAlET) p-FIUO*/> 
""'u 8*+3 M “ o»J 3 »b 


J-J3-M+1 

-tUTfc<to» 54 H* #I=i # 17 »a) 

foKMaTI l 2 tti* 2 > 
ir O^KiT.Eu.O) GO To 136 
/‘('ITctbtSAtil 

puh ^ r(/»2uy* , uu»j) values in a-adjacfnt fifl^*/) 

IuO-i 

L-L+l 

r up=±LoM5 

IML.Eg.6* IIJP=IL0+17 

"o I j2 *=1*J4 
J— J4*M+1 


’■h\1Tl<o»5CU> UH I# J) # 1=11.0* IUP) 
Fv.i 4 ' 1 ATTl ^PtJ* 1 ) 

’ w 0 =itlp 

.\ivlT4_(c*5^i> L 

^oK'IaT t//#10x* *N=* * !<J»//) 

Ir(L.EO.b> GO TO 12^ 

Gw To bOo 

EhL) 


ORIGINAL ? 

OF ROUE QcLiLal G 


68 


APPENDIX C 


TWO DIMENSIONAL ANALYSIS WITH 
DEPLETION APPROXIMATION 


< ^3 * 5 * ) 


310 


ttf 4 


*M"h*jIOn 
fuiCubiON 

•"iMENSION 
'*NiEivSiOM 
* * MEnS I ON 

HaMEnSION 9 ' 7 H r 9 uvn * vt » w ula 1 # 

UNFuSiON AC2{61)#BC^<bl> »0CZ(61 )#ECZl6i) 
n iMEuSiON a( 100) *CS<tot) #PHI (61) *ZflTt2) #GR<61) »*<93r6l » 
'-'i_ AL_ N 1 * N 1 2 
<^ADU6> h»V»U»UA 
i’i-wInU l ft 
'U=l.SttO 
fU2=NI + *2 


r'M4l* U) tU(S»’ 56) »UA( 1 3# 7) # 
aAR ( 26) »liAf<(26) *DAR<26) #EAR(26> 
AAZ12) »BA2U2) »UAZU2)»EAZ(l!>) 
AHr( 42) #Bl ! R<42> #UBR<4?) *EOR(42) 
AB2(12)fBr«2(l2) *UFiZU?> *EOZU?> 
ACR194) *BCR(94> *L»CK(94) *ECR(94) 


N«KIT=U 
Pt.Ar)( 5 » 804 ) 


Q|bt»QSS2 


format <2eio 

LuMAA-l 
L vhAA=l 
LuKAa=1 
Q=i .oE-59 
9bbl=QSSl*u 
0 bb2=0S52*u 
Vl=0.0259 
Ti-L=U*0 
ToA-0 • 4 
T J 1 =U • ti 
Cob=u • E 1 5 
tmuo- 2 • r^o 
•VTAl=3.9,c 
*i_TA<i-b«3c 
r>* = 3.l4l6 

rrSA=9*4*t.. 4b4E-i4 

Fh'bB- 1 1 • 7* 8 « 854F-14 

Et'bC=3» 7 *o * 4b4E“14 

FRbO=8*Bb4E-14 

U=9 

U=t/ 

TJ=?b 

Li=8 

l o — y 

K * — 1 U 

'•J * — 1 4 0 • 

1 4 --lb 2 * 

‘■0=140. 

L**=8 

14 -=^ 

t.J=32 
1 0=41 
Kt-=10 

Jc=U 

Lo=3* 

To=33 
I / = 9 * 
r/=93 
To=8*; 

'<o=5u 

.Jo=5i 

K*+=5b 

Kb=6u 

Jo — 6 X 

H=Tbl/( J3-3' 

CoMKjTl COEFFICIENTS In Nf-rlEU< 
r*v 2 1 = 1*11 
yU) = U-lJ^Wl/Ll 
nu 4 1=11 *1? 

Xa)=Wi + (l-Il)*W2/L? 
nu 6 1=12*13 
v l I ) + ( I-I 2 ) *W 3 /E 3 


jV ’- 1 

0 ? 


T5 


69 




in 


14 

16 


16 

?0 


2? 


C 


26 

pft 


20 

22 

2* 

36 


*j y. 


70 


Ju-T o~l 
PG ;i 1“?* f iU 

^mhm=i*7 

Pu lu J=l»Jl 

Vii)=-TSI*t ALPHA** )/<Al PHA-J*) 
nv, ur2#Kl 

a«z(v>>=fpwixij>-x< j-i> ) 

*MZiv()=rpbA»(XtJ+i}-xt^) )/d* 

<Wlm)=PAZ(J) , . 

-^Z(u3=FPbA/(XlJ+l)-x(J) ) 
rMLCO-MTE CO£FF iCIfcf ITS FOK V-UO-D. 
r, u m i— i * 14 

villi* i-1)*v:1/4 4.#L4)>3,*W1A* 

r u lo JL = I *+ r 7 . 
viI) = U-I**)*!rf?/L5*W1 
I o= I b”l 
f'o lev l=2»i'> 

fuK( il=Pl*tX(I+l)-X( 1-U3**XII + I )+2 f *X(T)+X(I-l> )/4, 
^(DsPIMHD+Xd-DF/fXdl-xa-ll) 

«-'ok( i)-PI*tv(I+l>+X(I> J/fxU + U-XllM 
r u M A >=ARMI) 

V w= J 

■iuPHh=1.2^ 

< aLPHA**K 2-1 *-K A* ( alpha- 1 . n* < alpha-i , 1 
*o= (ix2-l > +KLPHA*+K2-*2*ALPHA** 1^2-11 +1 • 

^rt=-AN/AO 

Ai_PHA=aLP«im+DA 
r-M] Xu=APSUjV ALPHA) 

Tr (PaVl^.OT.O.Ul) Tv '•U 

‘'u=T:»I/< 

-•ut jt=-HQM ALPHA* ♦( J2-J )-l*)/ *A< PHA-J# ) 

r 'w 24 J=2»h2 

AuZ(w>=FPbM/(?H|J)-7y{j-l) ) 

•uZ(o)=FPbH*(7ut J*-1)-Zb( J-l) )/2. 
r, uZ ( u 3 — Pbt 1 b ) 
ruZ l o) =FPbH/(?b l J+l) -Zb < f) ) 

C vt_Ft- 1 C I £) ■< 1 '■j FuP IJ— f iLLDi 
2u I“ 1 » 16 

/!!):( l-li+W2/(2.*L?*Lb)+*l 
‘ >u 7o I — Io # 1 7 

vi i )-XU6i + C 1- Ui)*wP/(2«*L5) 

Iu=I7“i 
'o 3u 

%«(l)=PI*(v{I + l)-X(I-l))*(XU + l )+2 i *XCl ) ♦X 4 1— 1 > )/ 4 . 
u.KU>=Pi*un>+x<i-i>)/um-xu-U > 

*'k-K ( 1 ) — PI * lX(I + l)+X(l) )/(XlI + H-X(D) 

ron C i ) =ACk ( T ) 

v k 1 ) =Zu ( J*--' 3 ) 

' k 2 ) -Zu < Jci-1 ) 
v k3)-0. 

’'u 2c 

vkJ)=TbT#tu-3)/tJ3-^> 

~'\j 34 J-o5 » J4 

v i J) -X l J3)+ToX*( J-J1)/(J4-J3> 

■’w 3u *J=rj4»Jb 

/ Vj)=XiJ4i+Ti:L*t J-J4 )/(J«5-J4> 

*cZl.J=FPW(Xl?)-XU) ) 

’’'-Z If) — FPbvt* t X 1 3) —X til) /? • 

rv,Z(4)-FPb*/tYl3)-XU> ) 

Av,z( j)-rpbM/(xi3)-x(c> ) • <F : c 

*'v4(^>-<FMi,P*(Xt4)-vi3) )+tPSA*{y(3)-x(2)) )/2. a- pnru. 

r ^vZ (o)-PCtV^) '-i rUOa 

f-\-z<o)-rpbo/(xi4)-x(o> > 

''o lo «j— 4 » k ^ 

^ wZ lo) =FPbu/ l X l^)) -X l J-l ) ) 

'cZl IX l J+1 )-a(J-1) )/2# 

Pv-4 l«j) -PC4i J) 
c uzto)-Fp^b/tyuu)-x(j) ) 

"oZ(tJ3)=FPbH/fXt J3)-X(J3-1) ) 

PL4lo3J=lt_HSC*lX(J3fl}-x(J3) >+LP5R*(X< J3)-X< J3-1 > ))/2# 
"CZ(J3)=HCZ< J3) 


^wzl j3>=EPbr/(Xl J3+' )-X< J3) > 
lurjH-1 


I 


i 

J 


*.ri 


>*? 


t. 3 

cion 
6i 2 




l;S 
i> b 

4» 

Ml 

F? 

M4 


r ;> l ie, O* luL f<Eb]ON. 

Ov 4 0 J=JL»JU 

J)=FRbC/(XJvf)-X( J-l) \ 

*H.*(u)=Ff>bt*<Xi J+l )->,( J> > /£!• 

<*C/(o)=PC<-U> 

roZ( J>=FPbL/(XU-H)-X(iJ> > 

'\vZ<o4)rtrS r /(X(J4)-A(J4-m 

<-4ZU4) = <i_F^ 0 *(y ( J4U )-X* J4) )+trsC*(Xt J4J-XC J4-1 ) ) >/2. 
'iwZtw4)=Bi-4( J4) 
fXZ(u4)=EHb r V(X( J4+1 )-X(J4>> 

.tu— 

■'u 4<^ J=JUrJU 
a^zc j)-rpbu/(xt J)-x(4-i) > 
it (u)=FPMJ*(XlJ+l)-Xld) )/2« 

' c/ u)-rpuo/(y(J+i)-/( j> ) 

C./'4''JTu 'Vjt'iMG A r i,> LIKETIWI. IHUT W 

' U H J vJ-~ ’j I ij 1 

M.i(u)=vi*bn/u*x ( j) 

-u( J)=u.h,+u!+rxf'{BLT£2*CT‘jl-XUi > i/TAUO 

• uA". VJ#VGf I>LRM»t WHIT 

r uN ,*T 3*2110) 

Tr ( ^j.U.U.n) GO TO blU 

’ ' u— >/v*“ A • 1 

’ iu i^Vu-2 . +VT* ALOG ( r f -> ( J 3 ) /UI ) 

AlU n O> Ji'iDA'-’Y VALUES* 

( 'W 44 1 = 1*11 

• i 1 » w 1 ) = V' J 
'u 4a 1=1 *a4 
’ ' M » vJ?> - Vw 

r w4ui — It - p T 3 
i 1 » oi J — n • 

"V 4C1 J ~ i * U * 

‘!llid)=Vj 
.u 17 » J)=pni ( J) 

u bij J=JJr J*4 

• i i .o)=Vj+ *(yr ~\iut /t J4-J3> 

.h 17# J)=Pnl ( J3) 

"u r 5t: 1 = 1 n R 
1 i i I * o4 ) =V b 

T^=l* IV-I4J 

"\j r >4 «j“l*b 

1 'H 1 1 t J ) =Vb 

,( TD* o)=rHT l J3) 


■ lu 'Au 


>« l = 

,b)=VG-(l-l)#(V -P«a< Ji>>/(T7-IA) 
I 1M tt-FirLD with KFLAXATIOM* 

I 


w>U . 

C r . * A < I 

*r mFG— u 

Tn<_o=l 

1 f G u CvhT*UUF 

rt'AGG=£PAbbU 

"..hT i — u 

T» ( T f j Ab^/ii^K.t.O,Lwt ) iT) IT-1 
EnPlH* ThK, GT/aTlMFuI CONTAIhr iMpLlLD hUiTIPl KATlCtl* 

Tv.i/j=<-1 )*♦< 1PASS) Q>-- b€.l±\*L 

t.^ti+i 

r oo2 I = T A » i n 
T>=4.*I~"**i 1 
i ? litul i-\l L ir>, J2> 

■■'••it-G b— 1 • H 

" i o~ g • rc : ' 1 - • - 1 ■* ” - ■■ • 

Lbbb-l Uj 


b4 


T uP=i3-1 
Jt_U=ui-l 
Jill =*i 

— l 

Tr ( IcUU.Eu* - 1 ) 

Jut)-L 
Jl I — 1 

llw— 1 

CbiiT i..iuF 


'w bci v»* Jl-Oi JHI.i 
o t»u 1=2 * 1' *P 
oL’l-WlIfOi 


0> 




H' 1 


GO 64 


71 


H 


7H 


5i ?. 


7 J? 
f 3( ! 6 
77 


7A 


(2 


■4 

f-5 


F.6 


*~aAiy< 1 ) ♦««.*»( J) 
i— bAtt 1 1 ) i J* 

"-uAn ( 1 ) *i ( J) 
r— uAivti ) *u«/t J) 

''-A+u+b+t 

i il-ia*w< i’J-i.i *3*u( i-i'j)+n*h< iM»j)+r*w< iiJ+iu/f 

\ l#u) = (i«-O l .iFbA)#y i ALD+oMfc'.( J A*V<TTt 

.= '/U* J)-.,ouO 


G > 
T'-' 




t>4 


‘ub^rtFS+O^ + S 
CoUTaHUF 
• ilf J>=WUf J) 

♦Ui3»J)=i«t IU>» J) 

r^NTiiiUF 
»'W 7o i = i t i3 
•iIfi)=wUfP) 
r wf4T xl IUF 
l uOP-Lo^P+1 
u (Loot .LL.L^PAX) 

Tr ( "^SiiT.t-U.O > GO 

.b=bOHT(KLF) 

(y! TtJ u » Fi/u ) 

r u K ! ««Tl/»£U Y » *UAJA r K ,i\ ,jU.J) KTFLDt,/) 
ulTLiu»5u<.l PM5 
r^Tl/.cUK ( *k^S=».| b.3»/> 
u 7 c T ■— 1 * 3 1 
.l-3l+l~ M 

i\ 1 Tl_ ( o » 4uu > ( 1 1 » J) * I-x» 
rubUNUF 
Fuk’^aTIIOa* 1 3F5. 

FwiiTiuur 

'V T.iE V-eirLP. 

i wLAi'E P0Ubn/ t P| VAUtS I* J AU; I* FTEtDF FMD Tr’+FRpOLAT^* 

T ,. = 1+3+11/4 

Tt^Icc-1 


>2) 


7o i = I«. * IK 
i/-4* t I-TW J +1 

V V K • i )=VLU 

anh »l> = wLO+<VHl«VL(j>/4* 
l ' lK + <;» 1 ) = VLO+ ( VHf-VLoi/2. 
tK + 3»l) = k/LFi+3. + ?VHT-VLO 

'/ IK + 4 » 1 ) = Wh! 

-ouTiUUF 

0^ J=2» J? 


VLOJ/4* 




F-17l>(G)-4o(1> >/7H(1 ) 


>-S*< V(Xl __ 

''UV#Jr=WU 0 #Kl>-S*<*tl2#Jl) 


f iMuF 

Tr ( IrAbF.ul ,1) GO TO 

»uO=14+P 
r :w Ah I = TLo# 1S 
'■-33+ A* ( T”iLO ) 

V l i>w(2)=UW<»j) 

CvuT i ! iUF 


« Jl)-w(lt »K1) ) . 


ttb 


V-FltLO, 


GO T' 1 flu 


■wNtiMUF 
' i_LAA l^t. 

< 'i*it_GA”i • £i 
I u UP - 1 
f uH- j.5-1 

tfc.0— u2“1 

Jt 1 1 -c. 

Jo=— i 

Ti t IbOi*. cj.-I 1 
Ji-O-c 

;-=l 
Cwl.Ti'iUF 

Pu 3— JL-U f 31 ! 1 * 3J 

'-■u Ao 

' * A Aj\ ( 1 ) +/ib? < J ) 
l i. 1 *bu7 ( J) 
til *u Li 7 ( J ) 

"-L4h( i)*ub7( J) 

O-A+d+u+F 

w vLO—V ( T t u t 

'J IL-(A*VU* J-D+fH-’tl-l* J>H)*VU*1»J)+F*V(I» J+1 >J /C 


ORlGHsAu 

0F Poor: ( J 




T 




■ i I iU> = < i .-OuFbA) *V« LLH0 M LGA*V1 T L 
'•w-V i 1 » J) ~w r *l_r' 
t-br.uL'jH)v**2 
?A CwIjt±;jof 
«>P CofjTifiUF 

»_oGP=LoOP+l 

rr ( luop.ll.lv^ax) to 46 
Tr l'\,RiT,t.u*Q> Go T n 95 
i\I Ti_ ( tj * 5Uo ) 

5ri8 rwh TaK/fcUY# *UATA FoR vU»J) FTFLD*/> 
i%lTL(b* ^Uc.) PMS 

Oc. H=1 *J? 

J-U2-M+1 

-<«IT^(b#5lU) (Vtlf J) ,I = i , I5#4) 

92 CuNTiNuF 
510 FvKMrtT ( 1 0 a r I lFd*2 ) 

P5 CwliTiNOF 

••'u ThE PtbAXAHOM Fop f>F U-FlEt D 

I.uAl) BOUrjuAPY VALUEF fro*’ THE V-FIEI.D ANu iNTF^POLATp 

1 1 1 1 t i, ) — V { 14 t JP “2 ) 

•111#*:)='/ {A4»J2“1) 

■h17»1)-VU5,J2-2) 

':vJ7,2)=VU5»Ji»-l) 

M117* 1)=VU4+1 »J2-2) 

*'u9o l=iO»T7»2 

‘*-l 4 * 2 +( I~i 6)/2 
Hilt i)=V(Hr Jl-2) 

Cl/nTINUF 
'•j 90 1 = 2*16 

! ‘ 4 1 * 1 > =' I U * I ) + 1 1 -1 ) * { iJ ( j 7 1 1 ) -U 1 1 « 1 M / 1 6 
ruUTiUlF 
”»w iuO t;xo»32 

i 1 1 1 * 1) =H( 1 7* 1) ♦ < 1-1^) # < u ( 33* 1 )»HC 17» 1) ) /if* 

roivlTINUF 
*uP=a7-1 

’*u 1 u2 T=J4» I UP *2 

'tlIul = (UU-l*l)tl|{ T+l, 1) >/2. 

ruliTiMLF 

!'Ub = U*U 

!.won=l 

0« - ic r m= 1 • H 
r uhT ^ilUF 

*♦*, ******«> 

CuMPIltu 

J-3 

702 

■ 4 lO — U • 

•'w i ~ u • 

* 1 1 —LI • 0 

Tt T'H I* J) .oF.UTEST) oO 700 

i*hO=ui*LF(o> 

J— uR t *• ) 

0 ”,«uL-l *F»o*ACPl I ) *H/2* 

T * J)=f)vO< *f r ’H04l.( r -4*AFR( I )*Ci r Sl 

" k i *«J >- n V^c + GFi^ 

V? '•iAFINIF 

w 7u6 J=h#K3 
r >j 7u6 I =«; r I UP 
' 4 i(J = 0 • 

,i_N=b * 

■ • I LSI =0.0 

T r I'M I * J) •G r .('TEST) oO TO 704 
■-’.lU = to *CSlU) ' 

’ l_l 1= IjK l B 1 

> 4 -« 0 L =1 .r-o*ACRin »h 

ul I » *J> =DXoL*RHO 
r ilio)=nvou»GFU 
fh ^RTA'iur 

j-j* 

r, w 7iU T=t.#TuP 

“I *U = u • 

' = U • 


96 


<-A 


U'O 


If 2 


Lf 4 



A 


’ 1 1 L ' * t = U • 0 

’t t U I #J) .G r .HTE5T) oO 70ft 
’i iU~\j*LS ( w J 
'-.uN=0U I J ) 

. » OL=L *F-o*ACP I 1 J ♦H/2« 

1 Jt I # J1 =L»VLrL*RHO + l , ('-4 *aCR < I ) *0 r S2 

"ijijl *=i'vwi.*cri4 


73 




^ vi '* T i i JUr 

r, t.LA a {Mr uU»J) FlrLD T l W006H FAPP, *GI1 . » AND rvjp t , 

" Aa = U • 

4.0= J 4-1 
Jf > i —c. 

1 

T r UcOO.Eu.-l) GO T« lub 

Jni =u4- 1 
Jw»— 1 

'v lib J — t JH1 t JJ 
Hu 114 i=<.#iur 
fl-AChU)*/iC7(J) 

-dC«U) *uCZ< J) 

'‘‘-iJCK J i > *uu? ( J j 
F-t-ChU >*«-(. 7( J> 
r-A+u+O+F 

■ I wLl- U t T # «j ) 

* iu=JA*u ,vj)+o^uf T+i.vJ) +r*u< i» J+u+o^d, j) )/r 

’ * i 1 § u ) - ( 1 • -Oj'iFOA ) +Uoi_D+o’4LbA*bl TL 
’r (MIT i J) .bT.rHI(J) ) u ( T » J ) “PHT ( J ) 

Ir ( J*Lb<* J^) GP TO 1^0 
Tu 11 0 

TrimlrJ) .lT.VIHV) > i ( I » j > =v inv 
1uLU-AdG<U( T r JJ-UOL”) 

Jr (AuUu.GI .Hu^-AX) O' lt.»AX=ALlUU 

CuWTit'lUF 

CwwTlNuF 

Th t^rtKiT.GT .U> WHlT r (b»bl4) DUMAX 
tl.o=ib +1 

Hu lib ! = ii.n f IUP 
^-mCk 1 1 ) *AC7( J4) 

P-UCht 1 ) *uL7( J4> 
n=uCMl)*ub7{j4) 

^-LCk 1 1 ) *t07< J4) 

C-A+u+U+E 

■ lyLO-Ut T r J41 
'*-1 + l-i A 

"i A*u t i» J.4-1 ) +b*ul i-1 t J4)+U*U< I + l* J4) +F*UA (m» 2 ) ) /C 
1 ‘ 1 1 » o4 ) =Yl »-0 M tGA ) *’ I0t-U+OMLGA*UTIL ‘ ' 

1)“UII »J4) 


AuU.J-AhS(UU r J4)-;iot ij) 

Jr (AuDU.Gi ,r>Uf"AX) U'VsAXsAObU 
rulMTiHUF 

1 ’>uLAa THE! uA—FIELO, 

Jt-O- J4+J 
Jt il=ub-l 

1^2 J=jlp»JHI 
l.-J-u4 + 1 

'u i t 0 T = iL r +IUP 
- 1 + 1 - J A 

^-hGk< J> 

1 1 ) *007 

M-UCKU )*ul 7( J) 
r -*cCn.t 1 1 *t»C? ( J) 
r-A+o+b+t 

' IwLHsUa ( m » l ) 

in lL- ( A*y^ir i .L*l)+b*uAlM-l»L)-»U*UA (^+l#L}+E*UA{-,L+l ))/r 
r. 1 N » L ) = ( 1 . -Of'tG A) ♦IJ01.D+0ME.GA*UTIL 
AuuMsAbS (UA< M#L)-U0LJ) 
ir (AoDU.Gl .^U'-AX) D’ iMAXsABUU 

G^MitiUF 

roKTlUUP 1 

Tr (.'*IUT,Ol ,0> WRITC'(6#5l^J DUMX 
r uH -SkT 1 4 OA t Fin * 3 ) 

r u+*0*T ( 2 Qa t r JO « 3) Oj?Tr*'?' r ^ 

LwuP^Lonp+i 

Tt (LuOP.Ll.LUMAX) &n jo '04 °I‘ > 

Tr L*,.RIT,G 1 .0) Go 5 J .7 m: ... 

■' 1 1 44 1 » J2J= U ( l 7»3) 

ruu=i4+? 

Ou U4 T-lLO,Ib 
-j7t2+( I-il.o) 

M Irw7J=UtN#bi 
Gv Tu mou 
rui^Ti'JuF 

rvi-iPuTL Ti it CHARGE r ,TOhl> IN I H r SIi I Com LAYFP. 

Cu -IrltTL Tl ir Cum^uT 
, ivp=u *5" 1 


CuK=0* 

OjIL= 0* 
ntLSui 

1 26 JS4* JUP 
nu 124 Ir2»IUP 
'iIL=QbTL + i.t-4*QS(I» J) 

'“wKsU.lh+Gt l*d> 

1 24 CuijTlNUF 
l?f* CwhTiUUF 

■'w i«i5 Ix2» TUP 

'•i-L=uEX*AU<< X)*tUt I » 1*K4> ) 
t?5 C^NTitJur 

Ou.L»l •L“4 : *‘t. r '!jC*GEU* ( J4«*J3) /TOX 
= i *u— 4*u+CLlK 

.i\ 1 Ti. t v # bio ) V J# VG » *5 IL » nO. » CUK 
516 Tuii‘.1nTt/»iU>*# *VJ=» #FlO,3,bX» 'VG=» 

♦ »*iLL=* #Fli,F,fjX**ISlL=* rElU.b) 

• *' ITu. 1 1> * 5 <jd ) 

5?6 JUV» * ABBREVIATE^ U-FIEIP*/| 

nw 5 m u M-1*J4*5 
J-J4-M 1 

540 ,kITc.(u»5hh) !si # I7#6) 

kITi_(o*53/> 

53? c’wK'-Otf l/* JOV, •AttiiREi'JArE’' OS-f*I r LD»/> 

5m2 ^=5#J3»4 

542 wnlTt<5#!>H4) <GS(It J) *I=1» |7#8> 

Wr^lTu t bt 5obl 

538 F wKMhT l /».iuv#* ABBREVIATED G-FIELD*/) 
2w 543 MSdfJ3f4 
JSJ3-M+3 

54 3 WKitt.(o#5H4) 

544 FuRMAti 1241 0.2) 

^KiTc.t5*54o) IFA5S 
546 FwHMATt/f-iuXF* !PAS5r**i4/l 
TuCPilNCPU 

iMlHAbS.tu.iTER^I G j To 800 
5u Tw moo 
61 0 CwNTiMUE 

*rvlTt^l6) w»V#U#UA 

OtftluD 16 

5»OP 

F,,u 


»Fin.3*bX»*0FT! =*,Flft,^,5V» 


ok-' ■"■ ,. 

QV 


"F U.S. GOVERNMENT PRINTING OFFICE 197W4V26E/393 REGION NO. 4 


75 


